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ABSTRACT 

We present initial results from the first systematic survey for Mg II quasar absorption lines at 
z > 2.5. Using infrared spectra of 46 high-redshift quasars, we discovered 111 Mg n systems over 
a path covering 1.9 < z < 6.3. Five systems have z > 5, with a maximum of z — 5.33 — the most 
distant Mg n system now known. The comoving Mg n line density for weaker systems (W r < 1.0A) 
is statistically consistent with no evolution from z = 0.4 to z = 5.5, while that for stronger systems 
increases three-fold until z ~ 3 before declining again towards higher redshifts. The equivalent width 
distribution, which fits an exponential, reflects this evolution by flattening as z — > 3 before steepening 
again. The rise and fall of the strong absorbers suggests a connection to the star formation rate density, 
as though they trace galactic outflows or other byproducts of star formation. The weaker systems' 
lack of evolution does not fit within this interpretation, but may be reproduced by extrapolating low 
redshift scaling relations between host galaxy luminosity and absorbing halo radius to earlier epochs. 
For the weak systems, luminosity-scaled models match the evolution better than similar models based 
on Mg II occupation of evolving CDM halo masses, which greatly underpredict dN/dz at early times 
unless the absorption efficiency of small haloes is significantly larger at early times. Taken together, 
these observations suggest that the general structure of Mg II -bearing haloes was put into place early 
in the process of galaxy assembly. Except for a transient appearance of stronger systems near the 
peak epoch of cosmic star formation, the basic properties of Mg II absorbers have evolved fairly little 
even as the (presumably) associated galaxy population grew substantially in stellar mass and half 
light radius. 

Subject headings: Galaxies: evolution — Galaxies: halos — Galaxies: high-redshift — Infrared: general — 
intcrgalactic medium — quasars: absorption lines 



1. INTRODUCTION 

Mg II absorption line searches provide a luminosity- 
independent method for studying the distribution 
of gas in galactic haloes. W hile this techniq u e has 
been employed for decades (IWevmann et al.l 119791: 
Lanzetta et al.l Il987t iTvtler et all 119871; iSargent et all 
19881: ISteidel fc Sargent! Il992t iNestor et all 120051: 
Prochter et al.l I2006D. us i ng t e ns of thousands o f 
sightlines (|Prochter etldl 120061: lOuider et all 120111 ). 
these systematic surveys have all focused exclusively 
on Mg II systems having z < 2.3. Only a handful of 
Mg II absorption systems have been reporte d at redshifts 
z > 3 based on serend i pitous detections ()Elston et all 
1996^ IKondo et al.ll200H I Jiang etalll200l . 

In this paper we present a systematic survey of inter- 
galactic Mg II absorption at 2 < z < 6 using QSO spectra 
taken with FIRE (the Folded-port InfraRed Echellette) 
on the Magellan Baade Telescope. By characterizing the 
evolution of Mg II absorption at higher redshifts we inves- 
tigate the mechanisms which populate extended haloes 
with magnesium during the epoch when the cosmic SFR 
approached and passed its peak at z « 2. 

Currently, two theories have been advocated to explain 
the properties of Mg II absorbers at low redshifts. In one 

1 This paper includes data gathered with the 6.5 meter Mag- 
ellan Telescopes located at Las Campanas Observatory, Chile 
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scenario, Mg II traces cool clumps embedded in heated 
outflows from galaxies with high specific star formation 
rates. Alternatively, the distribution of Mg n may re- 
flect gravitational and gas accretion processes, as from 
dynamical mergers or cold accretion filaments. 

The outflow hypothesis is supported by direct obser- 
vations of blueshifted Mg II absorption in t he spectra of 
star fo rming galaxies (jWeiner et all 120091 ): iRubin et all 
(2010) have also observed this trend and established a 
correlation between Mg n rest frame eq uivalent width 
and st ar formation rate (SFR). Moreover IBouche et all 
(2006) cross-correlated ~ 250,000 luminous red galax- 
ies with 1806 W r > 0.3A Mg 11 absorbing systems at 
z ~ 0.5 from SDSS-DR3, finding an anti-correlation be- 
tween Mg 11 rest frame equivalent width and galaxy halo 
mass. This suggests that the individual Mg 11 systems 
are not virialized, and the authors interpreted it as ev- 
idence that Mg 11 systems trace star formation driven 
winds. 

Concurrently, Zi betti et al.l (|2007f ) showed with a sam- 
ple of 2800 strong Mg 11 absorbers (W r > 0.8 A) at 
low redshifts (0.37 < z < 1.0) that W r and galaxy 
color are correlated, with stronger Mg II systems cor- 
responding to t he colors of blue star-forming galaxies . 
Follow up work (IBouche et alj|2007t [Menard et al.ll20llt 
iNoterdaeme et alJl2010HNestor et al.ll201lD supports this 
association of strong Mg II systems with star forming 
galaxies. 

Yet not all Mg n absorbers are found around star form- 
ing galaxies. In particular, galaxy-selected Mg II sam- 



2 



pies repeatedly show that star forming and early-type 
systems alike give rise to halo ab sorption, in propor - 
tions similar to the general field. iChen et al.l (|2010b) 
showed that the extent of the Mg n halo increases with 
galaxy stellar mass but correlates only weakly with spe- 
cific SFR for a sample of 47 weaker Mg n (mostly with 
W r < lA) at z < 0.5. The authors interpreted this as ev- 
idence that the Mg n absorbers reside in infalling clouds 
that eventually fuel star formation. In contrast to Zi- 
betti, Ch en et al.1 (|2010a[ ) see little correlation between 
W r strength and galaxy colors, a r esult corroborated by 
later research on weak abs o rbers (lLovegrove fc Simcod 
[20Tlt iKacprzak et al.|[20TTI ). iKacprzak et all (|2011[ ) re^ 
port a correlation between W r and galaxy inclination in 
a sample of 40 Mg II absorption selected galaxies with 
0.3 < z < 1 and W r < lA. This hints that some 
Mg II absorbers may exhibit co-planar geometries and 
organized angular momenta, as may be found in accret- 
ing streams, satellites, and filaments, in contrast to star 
formation driven winds, which escape perpendicular to 
the disk. 

The different interpretations of the above studies may 
be understood in terms of the selection processes used 
to construct each sample. Broadly speaking, galaxy 
surveys around QSO sightlines pre-selected for strong 
Mg II absorption tend to find systems with high specific 
star fo rmation rates dZibetti et al.1 120071 iBouche et all 
20071 iMenard et all 120111: iNoterdaeme et al.l 120101: 
Nestor et al.l 120 111) . Conversely, blind searches for 
Mg II absorption along QSO sightlines serendipitously 
located near galaxies tend to find weaker Mg n systems, 
and little correlatio n between Mg II incidence a nd SFR 
(jChen et alll2010bl lal: lLovegrove fc SimcodfeOllI ). 

Recently, the idea of a bimodal Mg n absorber 
population has gained favor, with weaker absorbers 
W r < lA mostly tracing distributed halo gas and 
stronger absorbers W r > l A mostly tracing outflows 
(Kacprza k fc Churchilll 120111 ). This naturally explains 
the discrepant results between Mg n-selected galaxies 
and galaxy-selected Mg II absorbers: the two meth- 
ods simply sample different ends of the Mg n equiva- 
lent width distribution. While this explanation is sim- 
ple and appealing, the existence of strong absorbers 
around fairly quiescent galaxies dGauthier et all 120101 : 
IBowen fc Chelouchd l20Tlt IGauthier fc Chenll201lD sug- 
gests that in reality there may be significant overlap be- 
tween the populations in W r space. 

The volume-averaged SFR, halo assembly and cold ac- 
cretion rates, and metal enrichment rate all vary with 
lookback time, but the restriction of z < 2 for optical 
Mg II searches means that these surveys do not sample 
epochs predating the SFR peak at z ~ 2.5 — 3. This in- 
flection in the SFR density, along with the shutdown of 
cold accretion at low redshift, could be a useful diagnostic 
for evaluating the plausibility of the various mechanisms 
for distributing Mg II into haloes. 

While these factors motivate the study of Mg II systems 
at high redshift, the endeavor has proven difficult in prac- 
tice since the observed-frame transition moves into the 
near-infrared. In this regime, absorption line searches are 
hampered by blending with OH sky emission and telluric 
absorption, except at high resolutions where instrument 
sensitivity becomes an issue. We recently commissioned 



the FIRE infrared spectrometer on Magellan with the 
aim of studying C iv at z > 6; its design characteristics 
also make it a powerful survey instrument for mapping 
Mg II between z ~ 2 and z ~ 5.5 (and beyond as back- 
ground targets are discovered). This paper describes the 
first results of an ongoing study of Mg II absorption in 
early universe, based on an initial sample of 46 sightlines 
observed with FIRE. 

In Section [2] we describe the data acquisition. In Sec- 
tion [3] we describe our analysis, including our Mg n line 
identification procedure and completeness tests. In Sec- 
tion [4] we provide our main science results, including 
binned and maximum likelihood fit population distri- 
butions, and compare these results to those of previous 
studies with lower redshift search ranges. In Section [5j 
we discuss the implications of these results for the origin 
of Mg II absorbing systems. Throughout this paper we 
use a ACDM cosmology with ft m — 0.3, — 0.7, and 
H Q = 70 km s" 1 Mpc-\ 

2. DATA ACQUISITION 

Our sample consists of 46 qu asar spectra taken with 
FIRE (jSimcoe et al.ll2008ll20Toh . between 2010 June and 
2011 April. This instrument is a single object, prism 
cross-dispersed infrared spectrometer with a spectral res- 
olution of R=6000, or approximately 50 km/s, over the 
range 0.8 to 2.5 /im. We observed using a 0.6" slit and 
with a typical seeing of 0.5"-0.8" . The wavelength cover- 
age of FIRE set the lower redshift limit for our Mg n ab- 
sorption doublet search at z « 1.9 in each sightline. The 
QSO emission redshifts, which ranged from z — 3.55 
to z — 6.28, typically set the upper limits (although 
masked telluric regions described below sometimes de- 
creased these maximum cutoffs). Regions lying within 
3000 km s" 1 of each background QSO were excluded 
from the search path, although the effect of this cut was 
negligible on our derived results. Table [T] lists key obser- 
vational properties of the 46 quasars. The quasars them- 
selves were pred ominantly chosen from the SDSS DR7 
quasar catalog (jSchneider et al.l l2010f ) . with an addi- 
tional, somewhat heterogeneous selection of bright, well 
known objects from the literature. 

We reduced the data using a custom-developed IDL 
pipeline named FIREHOSE, which has been released to 
the community as part of the FIRE instrument pack- 
age. This pipeline evolved from the optical echelle re - 
duction software package MASE ( Bochan ski et al.l l2009) . 
FIREHOSE contains many important modifications for 
IR spectroscopy, including using terrestrial OH lines im- 
printed on science spectra for wavelength calibration. 

To correct for telluric absorption features, we obtained 
spectra of A0V stars at comparable observing times, air- 
masses and sky positions to our observed QSOs. We then 
performed telluric corrections and relative flux calibra- 
tions using the xtellco r software package r eleased with 
the spextool pipelin e dCushing et al.l I2004D . which em- 
ploys the method of iVacca et all (|2003f ). Regions lying 
between the J and H bands (1.35 /Ltm to 1.48 fim, ex- 
cluding z =3.8 to 4.3 for Mg II systems) and the H and 
K bands (1.79 fim to 1.93 ^m, excluding z =5.4 to 5.9), 
were masked out of the search path. Another smaller 
region of path length Az ~ 0.2 centered on z ~ 3.1 
(near the Y-J transition) showed poor telluric residu- 
als (noticeable in Figure Q} . We left this region in our 
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Figure 1. Three example spectra arranged from top to bottom in order of decreasing signal- to- noise ratio (SNR). The top spectrum is 
the highest-quality in the sample (SNR= 47), the middle is at the survey median quality (SNRfa 13) and the bottom is one of the poorest 
quality spectra (SNR= 5). The Mg II emission line is clearly visible at A reat ss 2800A, as is a noisy region from poor telluric correction 
centered at ~ 1.14 fim in the observed frame (z ~ 3.1 for Mg II absorbing systems.) 

our highest (top panel), a typical (middle panel), and a 
relatively poor (bottom panel) signal-to-noise spectra. 
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Figure 2. Representative examples of Mg II absorption from the 
survey. These are chosen to illustrate some of the common features 
in the infrared data, including blending with sky emission lines, and 
variable signal-to-noise ratios. Several very strong, and saturated 
Mg II systems are present in the data, sometimes with complex 
velocity structure. In the bottom right panel, we see a pair of 
Mg II doublets with blended 2796/2803 components, separated by 
750 km/s. 



Mg II search since it was still possible to find strong sys- 
tems (we found six) . We report these systems in our line 
lists but exclude this path and associated systems from 
our population statistics as the sample completeness is 
low and variable from one sightline to the next. 

The final signal-to-noise ratios per pixel for our full 
QSO set ranged from 4.0 to 47.2, with a median value 
of 12.9. Figure Q] shows the Mg n search portions for 



3. ANALYSIS 
3.1. Mg II Line Identification 
3.1.1. Continuum Fitting 

With just 46 QSO sightlines in our sample, it would 
be reasonable to fit each continuum manually. However, 
as part of our downstream analysis we also run extensive 
Monte Carlo tests to characterize sample completeness. 
We wished to capture the effects of the continuum fitting 
on our completeness. Accordingly, we developed an au- 
tomated algorithm for continuum fitting, to make treat- 
ment of the Monte Carlo sample both tractable, and as 
consistent as possible with our handling of the true data. 

First, we generated an initial line mask to prevent ab- 
sorption and emission features (including miscorrected 
skylines) from biasing the continuum. Next, we per- 
formed iterative, sigma-clipped linear fits on the masked 
data over small segments of the spectra (~100 pixels, or 
^1250 km/s, wide) in order to determine a list of spline 
knots (with two knots drawn from each of these small seg- 
ments). We then implemented a cubic spline fit through 
these knots to produce a continuum fit. Since knowledge 
of the continuum fit facilitates the creation of the initial 
line mask, we iterated this between two and five times 
before converging on a final continuum fit. 

3.1.2. Matched Filter Search 
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Table 1 

FIRE Mg II Survey Sightlincs 



Object 


z QSO 


Az 


texp 

« 


SNR a 


Q0000-26 


4.10 


1.95-3.83 


1226 


20.7 


BR0004-6224 


4.51 


1.95-4.51 


1764 


8.2 


BR0016-3544 


4.15 


1.95-3.83 


2409 


14.0 


SDSS0106+0048 


4.45 


1.95-4.45 


3635 


18.9 


SDSS0113-0935 


3.67 


1.95-3.67 


1944 


12.8 


SDSS0127-0045 


4.08 


1.95-3.83 


3635 


22.5 


SDSS0140-0839 


3.71 


1.95-3.71 


1226 


18.2 


SDSS0203+0012 


5.85 


1.98-5.40 


3635 


4.0 


SDSS0244-0816 


4.07 


1.95-3.83 


1944 


12.9 


BR0305-4957 


4.78 


1.95-4.78 


2409 


28.2 


BR0322-2928 


4.62 


1.95-4.62 


2409 


21.1 


BR0331-1622 


4.32 


1.95-4.32 


1944 


15.1 


SDSS0331-0741 


4.74 


1.95-4.74 


2177 


6.2 


SDSS0332-0654 


3.69 


1.95-3.69 


2409 


5.6 


SDSS0344-0653 


3.96 


1.95-3.83 


3022 


6.6 


BR0353-3820 


4.58 


1.95-4.58 


1200 


26.7 


BR0418-5723 


4.37 


1.95-4.37 


4200 


8.5 


SDSS0818+1722 


5.90 


2.00-5.40 


9000 


10.2 


SDSS0836+0054 


5.82 


1.96-5.40 


10187 


15.8 


SDSS0842+0637 


3.66 


1.95-3.66 


2409 


9.1 


SDSS0935+0022 


5.82 


1.96-5.40 


1817 


11.1 


SDSS0949+0335 


4.05 


1.95-3.83 


1817 


13.6 


SDSS1020+0922 


3.64 


1.95-3.64 


2409 


15.2 


SDSS1030+0524 


6.28 


2.16-6.28 


14400 


5.0 


SDSS1037+0704 


4.10 


1.95-3.83 


2726 


8.1 


SDSS1110+0244 


4.12 


1.95-3.83 


2409 


18.6 


SDSS1135+0842 


3.83 


1.95-3.83 


2409 


17.7 


SDSS1249-0159 


3.64 


1.95-3.64 


1817 


18.2 


SDSS1305+0521 


4.09 


1.95-3.83 


1363 


8.8 


SDSS1306+0356 


5.99 


2.04-5.99 


15682 


6.4 


ULAS1319+0950 


6.13 


2.10-6.13 


19275 


5.0 


SDSS1402+0146 


4.16 


1.95-3.83 


1902 


15.0 


SDSS1408+0205 


4.01 


1.95-3.83 


2409 


9.9 


SDSS1411+1217 


5.93 


2.01-5.93 


3600 


8.6 


01 422-4-2309 


3.65 


1 95-3 65 


1226 


47.2 


SDSS1433+0227 


4.72 


1.95-4.72 


2409 


13^4 


CFQS1509-1749 


6.12 


2.10-6.12 


9900 


9.7 


SDSS1538+0855 


3.55 


1.95-3.55 


1363 


24.2 


SDSS1616+0501 


4.87 


1.95-4.87 


1800 


5.2 


SDSS1620+0020 


4.09 


1.95-3.83 


972 


7.0 


SDSS1621-0042 


3.70 


1.95-3.70 


1204 


26.1 


SDSS2147-0838 


4.59 


1.95-4.59 


2409 


13.8 


SDSS2225-0014 


4.89 


1.95-4.89 


7200 


11.5 


SDSS2228-0757 


5.14 


1.95-5.14 


3600 


4.9 


SDSS2310+1855 


6.04 


2.06-6.04 


14400 


17.5 


BR2346-3729 


4.21 


1.95-3.83 


2409 


11.0 


a Median signal-to-noise 
Mg II pathlength. 


ratio per 


pixel 


across 



After obtaining continuum-normalized spectra, we ran 
a matched filter search with a signal-to-noise ratio cut 
of 5 to determine an initial candidate list. The filter 
kernel consists of a pair of Gaussian profiles separated 
by the Mg 11 spacing with FWHM ranging from 37.5 
to 150 km/s. Sky subtraction residuals near OH lines 
are a source a noise in the infrared, so we identified and 
masked badly miscorrected sky lines before running this 
convolution to reduce the number of false positives. Even 
with this masking, this matched filter search typically 
led to ^50 to 100 candidates per QSO sightline, of which 
only ~ 2 per sightline were true detections. 

To mitigate the high false positive rate of the matched 
filter procedure, we fit each candidate doublet with a pair 
of Gaussian absorption profiles and subjected the result- 
ing fit parameters to a set of consistency checks (e.g. 
W2796 > W2803 within errors). This additional doublet 
fitting filter step reduced the number of false positives by 
a factor of ~10, and left us with 256 candidates across 



all QSO sightlines. 

3.1.3. Visual Inspection 

We subjected these machine generated candidates to 
a visual inspection to produce the final Mg 11 sample 
list. This interactive step lowered the final list to 110 
Mg 11 systems (plus one proximate system) . Despite our 
best efforts to automate the process, we found that the 
poorly telluric corrected regions and abundance of sky 
lines mandated this visual inspection. 

We also visually combed through all 46 spectra to 
identify systems missed by the automated line finder. 
We found three systems with this visual search that our 
automated line finding algorithm rejected. One system 
was a large 2.62A system that was "too complex" (and 
therefore was fit poorly by a Gaussian absorption pro- 
file and consequently rejected), and the other two were 
smaller systems (confirmed by Fell lines) that were par- 
tially obscured by sky lines. The two smaller systems 
were manually added to our line list in Table [21 but were 
not included in statistical calculations so as not to bias 
our results by counting them and then overcorrecting for 
completeness. The large system's equivalent width was 
greater than the highest value (W r = 0.95A) probed in 
our Monte Carlo completeness tests. We are ~ 100% 
complete to such large systems, so this system was in- 
cluded in our statistical calculations. 

3.1.4. Final Mg 11 Sample 

The final Mg 11 line list, as well as the main properties 
of each system, are listed in Table [2] Figure [2] shows a 
selection of six Mg 11 systems in detail, while Figure [3] 
shows continuum normalized postage stamps of the en- 
tire sample. 

3.2. Simulations of Completeness and Contamination 

The nature of the infrared sky suggests that our sur- 
vey completeness will vary strongly with redshift as sys- 
tems overlap with atmospheric OH features. The com- 
mon method of calculating sensitivity from ID error vec- 
tors, usually estimated assuming photon statistics of the 
background, may not capture sky subtraction residuals 
completely. For this reason we have performed extensive 
simulations to characterize completeness and contamina- 
tion in both the automated and interactive portions of 
our line identification procedure. 

3.2.1. Automated Completeness Tests 

To test the automated portion of our identification 
pipeline, we ran a large Monte Carlo simulation in which 
we injected artificial Mg 11 systems into our QSO spec- 
tra and then ran our line identification algorithm to de- 
termine success/recovery rates. The full Monte Carlo 
completeness simulation involved 20,000 simulated spec- 
tra per QSO injected with Mg 11 at an inflated rate of 
dN/dz « 5, which combined led to over 12 million in- 
jected systems. The velocity spreads of the injected sys- 
tems were calibrated to mimic the behavior of true, de- 
tected systems (see Figure [Til below) . The injected sys- 
tems were distributed uniformly in rest frame equivalent 
between 0.05A < W < 0.95A. (A separate test showed 
that the completeness did not improve significantly be- 
yond this upper limit.) Our small number of QSOs al- 
lowed us to calculate individual completeness matrices 
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Figure 3. The full FIRE Mg II sample, in continuum-normalized units, ovcrplotted with errors (dotted line). The bottom edge of each 
plot corresponds to 0, and the horizontal line near the top to 1. Panel numbers correspond to entries in Table [2] 



6 



Table 2 

Summary of Absorption Properties for the FIRE Mg II Sample 
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i 


Q0000-26 


2.184 


0.162 


0.031 


0.112 


0.026 


166.0 


2 


Q0000-26 


3.390 


1.356 


0.016 


1.145 


0.016 


259.7 


3 


BR0004-6224 


2.663 


0.260 


0.045 


0.140 


0.036 


58.0 


4 


BR0004-6224 


2.908 


0.596 


0.047 


0.183 


0.028 


83.3 


5 


BR0004-6224 


2.959 


0.569 


0.063 


0.669 


0.047 


108.4 


6 


BR0004-6224 


3.203 


0.558 


0.029 


0.548 


0.026 


86.8 


7 


BR0004-6224 


3.694 


0.236 


0.042 


0.234 


0.019 


79.6 


8 


BR0004-6224 


3.776 


1.045 


0.046 


1.049 


0.043 


169.7 


9 


BR0016-3544 


2.783 


0.517 


0.027 


0.305 


0.026 


101.2 


10 


BR0016-3544 


2.819 


4.028 


0.050 


3.639 


0.053 


605.5 


11 


BR0016-3544 


2.949 


0.157 


0.026 


0.144 


0.035 


79.6 


12 


BR0016-3544 


3.757 


1.559 


0.041 


1.430 


0.050 


551.4 


13 


SDSS0106+0048 


3.729 


0.842 


0.016 


0.673 


0.015 


194.9 


14 


SDSS0113-0935 


2.825 


0.194 


0.029 


0.110 


0.029 


94.0 


15 


SDSS0113-0935 


3.544 


0.228 


0.037 


0.186 


0.035 


133.6 


16 


SDSS0113-0935 


3.617 


0.563 


0.024 


0.344 


0.020 


155.3 


17 


SDSS0127-0045 


2.588 


1.602 


0.025 


1.164 


0.025 


457.8 


18 


SDSS0127-0045 


2.945 


2.253 


0.038 


1.583 


0.037 


529.7 


19 


SDSS0127-0045 


3.168 


0.309 


0.024 


0.138 


0.016 


198.5 


20 


SDSS0127-0045 


3.728 


0.824 


0.012 


0.745 


0.012 


176.9 


21 


SDSS0140-0839 


2.241 


0.405 


0.031 


0.686 


0.043 


144.4 


22 a 


SDSS0140-0839 


3.081 


0.558 


0.018 


0.410 


0.027 


126.4 


23 


SDSS0140-0839 


3.212 


0.081 


0.014 


0.092 


0.014 


126.4 


24 b 


SDSS0140-0839 


3.697 


0.604 


0.014 


0.308 


0.009 


173.3 


25 


SDSS0203+0012 


3.711 


0.374 


0.038 


0.250 


0.065 


122.8 


26 a 


SDSS0203+0012 


4.313 


0.849 


0.093 


0.824 


0.080 


198.5 


27 


SDSS0203+0012 


4.482 


0.670 


0.183 


0.623 


0.024 


126.5 


28 


SDSS0203+0012 


4.978 


0.886 


0.039 


0.791 


0.056 


151.6 


29 


BR0305-4957 


2.502 


0.331 


0.024 


0.169 


0.022 


137.2 


30 


BR0305-4957 


2.629 


1.113 


0.018 


0.959 


0.023 


220.1 


31 


BR0305-4957 


3.354 


0.564 


0.013 


0.412 


0.012 


144.5 


32 


BR0305-4957 


3.591 


1.373 


0.017 


1.207 


0.013 


306.5 


33 


BR0305-4957 


4.466 


1.792 


0.017 


1.478 


0.029 


302.9 


34 


BR0322-2928 


2.229 


0.618 


0.020 


0.510 


0.021 


112.1 


35 


BR0331-1622 


2.295 


1.836 


0.067 


1.714 


0.056 


403.8 


36 


BR0331-1622 


2.593 


0.223 


0.019 


0.185 


0.019 


76.0 


37 


BR0331-1622 


2.927 


1.382 


0.039 


1.098 


0.045 


346.1 


38 


BR0331-1622 


3.557 


0.707 


0.033 


0.582 


0.033 


130.0 


39 a 


SDSS0332-0654 


3.061 


0.883 


0.084 


0.608 


0.059 


162.4 


40 


BR0353-3820 


1.987 


3.131 


0.030 


2.717 


0.026 


508.2 


41 


BR0353-3820 


2.696 


0.381 


0.014 


0.232 


0.014 


187.7 


42 


BR0353-3820 


2.754 


4.599 


0.016 


4.325 


0.019 


616.3 


43 


BR0418-5723 


2.030 


1.449 


0.072 


1.009 


0.080 


209.3 


44 a 


BR0418-5723 


2.978 


1.850 


0.072 


2.136 


0.099 


263.3 


45 


SDSS0818+1722 


3.563 


0.640 


0.072 


0.427 


0.029 


151.7 


46 


SDSS0818+1722 


4.431 


0.457 


0.052 


0.138 


0.010 


108.4 


47 


SDSS0818+1722 


5.065 


0.841 


0.061 


0.533 


0.046 


151.7 


48 


SDSS0836+0054 


2.299 


0.455 


0.022 


0.300 


0.021 


130.0 


49 


SDSS0836+0054 


3.744 


2.607 


0.024 


1.992 


0.031 


443.4 


50 


SDSS0949+0335 


2.289 


2.852 


0.062 


2.408 


0.054 


421.8 


51 


SDSS0949+0335 


3.310 


2.033 


0.039 


1.665 


0.033 


339.0 


52 


SDSS1020+0922 


2.046 


0.106 


0.045 


0.288 


0.047 


112.0 


53 


SDSS1020+0922 


2.593 


0.464 


0.026 


0.499 


0.022 


144.5 


54 


SDSS1020+0922 


2.749 


0.652 


0.023 


0.518 


0.024 


151.7 


55 


SDSS1020+0922 


3.479 


0.118 


0.016 


0.085 


0.019 


104.8 


56 


SDSS1030+0524 


2.188 


0.317 


0.017 


0.291 


0.017 


133.7 


57 c 


SDSS1030+0524 


2.780 


2.617 


0.069 


1.855 


0.086 


583.9 


58 


SDSS1030+0524 


4.583 


1.857 


0.031 


2.139 


0.127 


367.7 


59 


SDSS1030+0524 


4.948 


0.447 


0.017 


0.278 


0.019 


162.4 


60 


SDSS1030+0524 


5.130 


0.138 


0.013 


0.089 


0.023 


76.0 
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Table 2 

Summary of Absorption Properties for the FIRE Mg II Sample (Continued) 



Index # Sightline z W/ r (2796) <r(2796) W r (2803) cr(2803) Av 











(k) 
( A ) 




(k\ 
\ A ) 


\ A ) 


^Km s J 


61 


SDSS1110H 


-0244 


2.119 


3.041 


0.041 


2.884 


0.042 


454.2 


62 


SDSS1110H 


-0244 


2.223 


0.205 


0.024 


0.121 


0.029 


166.0 


63 


SDSS1305H 


-0521 


2.302 


1.993 


0.095 


1.533 


0.095 


234.5 


64 


SDSS1305H 


-0521 


2.753 


0.378 


0.040 


0.319 


0.037 


140.9 


65 


SDSS1305H 


-0521 


3.235 


0.328 


0.025 


0.130 


0.025 


122.8 


66 


SDSS1305H 


-0521 


3.680 


1.781 


0.068 


1.583 


0.039 


353.4 


67 


SDSS1306H 


-0356 


2.533 


3.307 


0.101 


3.019 


0.088 


594.6 


68 


SDSS1306H 


-0356 


3.490 


0.648 


0.031 


0.526 


0.044 


169.7 


69 


SDSS1306H 


-0356 


4.615 


0.983 


0.078 


0.724 


0.038 


187.7 


70 


SDSS1306H 


-0356 


4.865 


2.798 


0.044 


3.049 


0.087 


403.7 


71 


SDSS1306H 


-0356 


4.882 


1.941 


0.079 


2.276 


0.040 


248.8 


72 


ULAS1319 


4-0950 


4.569 


0.406 


0.062 


0.177 


0.029 


101.2 


73 


SDSS1402H 


-0146 


3.277 


1.075 


0.018 


1.034 


0.028 


234.4 


74 c 


SDSS1402+0146 


3.454 


0.341 


0.016 


0.112 


0.018 


173.3 


75 


SDSS1408H 


-0205 


1.982 


2.174 


0.056 


1.769 


0.054 


252.5 


76 


SDSS1408H 


-0205 


1.991 


0.914 


0.041 


0.555 


0.042 


126.4 


77 


SDSS1408H 


-0205 


2.462 


1.385 


0.040 


1.029 


0.035 


209.3 


78 


SDSS141H 


-1217 


2.237 


0.627 


0.041 


0.334 


0.045 


148.0 


79 


SDSS141H 


-1217 


3.477 


0.343 


0.016 


0.179 


0.022 


83.2 


80 


SDSS141H 


-1217 


4.929 


0.644 


0.023 


0.488 


0.018 


317.2 


81 


SDSS141H 


-1217 


5.055 


0.207 


0.013 


0.092 


0.015 


90.4 


82 


SDSS141H 


-1217 


5.250 


0.330 


0.013 


0.190 


0.011 


130.1 


83 


SDSS1411H 


-1217 


5.332 


0.197 


0.013 


0.241 


0.011 


97.6 


84 c 


Q1422+2309 


3.540 


0.169 


0.011 


0.115 


0.008 


130.0 


85 


SDSS1433H 


-0227 


2.772 


0.735 


0.018 


0.601 


0.024 


166.0 


86 a 


CFQS1509-1749 


3.128 


0.858 


0.093 


0.773 


0.043 


220.1 


87 


CFQS1509-1749 


3.266 


0.896 


0.021 


0.711 


0.023 


194.9 


88 


CFQS1509-1749 


3.392 


5.585 


0.071 


5.082 


0.050 


709.8 


89 


SDSS1538H 


-0855 


2.638 


0.278 


0.027 


0.206 


0.028 


180.4 


90 


SDSS1538+0855 


3.498 


0.151 


0.011 


0.122 


0.014 


173.2 


91 


SDSS1616H 


-0501 


2.741 


1.510 


0.044 


0.923 


0.051 


169.7 


92 


SDSS1616H 


-0501 


3.275 


0.600 


0.036 


0.494 


0.110 


104.8 


93 


SDSS1616H 


-0501 


3.396 


0.960 


0.036 


0.631 


0.113 


126.4 


91 


SDSS1616H 


-0501 


3.450 


0.606 


0.033 


0.557 


0.053 


115.6 


95 


SDSS1616H 


-0501 


3.733 


2.252 


0.189 


1.421 


0.068 


263.3 


96 


SDSS1620H 


-0020 


2.910 


1.130 


0.058 


1.063 


0.058 


263.3 


97 


SDSS1620H 


-0020 


3.273 


0.965 


0.043 


0.635 


0.052 


202.1 


98 


SDSS1620H 


-0020 


3.620 


1.357 


0.065 


1.091 


0.042 


320.9 


99 


SDSS1620H 


-0020 


3.752 


1.656 


0.065 


1.550 


0.095 


374.9 


100 


SDSS1621-0042 


2.678 


0.176 


0.017 


0.135 


0.016 


97.6 


101 a 


SDSS1621-0042 


3.106 


0.974 


0.011 


1.011 


0.012 


112.0 


102 


SDSS2147-0838 


2.286 


0.977 


0.040 


0.567 


0.033 


158.9 


103 


SDSS2228-0757 


3.175 


0.304 


0.037 


0.243 


0.031 


50.8 


104 


SDSS2310+1855 


2.243 


1.441 


0.050 


0.781 


0.049 


292.1 


105 


SDSS2310+1855 


2.351 


0.807 


0.042 


0.492 


0.035 


212.9 


106 


SDSS2310+1855 


2.643 


0.863 


0.036 


0.339 


0.059 


173.3 


107 


SDSS2310+1855 


3.300 


0.665 


0.039 


0.457 


0.034 


212.9 


108 


BR2346-3729 


2.830 


1.633 


0.049 


1.421 


0.037 


238.1 


109 


BR2346-3729 


2.923 


0.557 


0.030 


0.636 


0.034 


162.5 


110 


BR2346-3729 


3.619 


0.412 


0.031 


0.240 


0.019 


151.6 


111 


BR2346-3729 


3.692 


0.385 


0.016 


0.413 


0.046 


133.6 



a Poor telluric region 
b Proximate system 

c Missed by automated search algorithm 

for each sightline q, in finely tuned redshift (dz = 0.02) 
and rest frame equivalent width (dW — O.OlA) bins. We 
denote these values, which reflect the completeness of the 
automated line search alone, as L q (z, W). 

3.2.2. False Positive and Spurious Rejection Tests 

Because our identification algorithm involves an inter- 
active evaluation of each system, the user may reject 
true Mg II systems or accept spurious doublet candidates 
caused by miscorrected skylines or correlated noise. We 
developed a simple test to quantify these errors using 
simulated data. First, we created a large number of sim- 
ulated spectra in which we injected Mg n doublets us- 
ing a process identical to that described in Section 13.2.1 1 



above, but with a slightly larger doublet velocity spac- 
ing. We ran these artificial spectra through our identifi- 
cation and sample cleaning algorithms, and interactively 
accepted/rejected over 1500 such candidates. By using 
nonphysically spaced doublets, we were guaranteed that 
auto-identified candidates were either correlated noise or 
simulated Mg n doublets. The human decision accuracy 
for each of these cases was translated into success rates, 
which were then folded into our completeness results. 

The time consuming nature of this interactive test did 
not allow for the fine grained calculations produced above 
in the automated completeness tests. We found, how- 
ever, that our ability to distinguish both false positives 
and true systems scaled strictly with the signal-to-noise 
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Figure 5. The pathlcngth-weighted average completeness 
C(z, W) as a function of redshift z and equivalent width W r , for all 
survey sightlines combined. The error in this estimate is typically 
< 5% within our the search range. 



Figure 4. Calibration of user accuracy for interactive acceptance 
of Mg II doublets and rejections of false-positive candidates. The 
individual points represent the binned rates, and the overplottcd 
solid lines are the maximum likelihood fits to the functional forms 
in Equations [T] (top panel) and [2] (bottom panel). Displayed errors 
bars on the binned acceptance rates represent the Wilson score 
interval. 

ratio of the detected candidates. We therefore used this 
simulation to parametrize these success rates as a func- 
tion of line SNR, and applied the resulting scalings to 
the larger Monte Carlo completeness simulation and true 
data sets. 

Specifically, for each candidate we calculate a boxcar 
SNR s = W r /a Wr of the 2796 A line. We then tab- 
ulated the user-accept percentages of both the injected 
systems (P MgII (s); ideally 100%), and the false positives 
(P FP (s), ideally 0%) as a function of the SNR by cal- 
culating maximum likelihood fits to chosen functional 
forms. For the injected systems, we chose an exponential 

for pMgii 



P Mg u (s) = P^ (l - e -''°) , 



(1) 



where Px, (MLE estimate: 0.984) is the probability that 
the user accepts a true system with large SNR and s c 
(MLE estimate: 2.89) is an exponential scale factor. 
Note that even at large SNR the acceptance rate is not 
100% because of residual regions with very bright sky 
lines and/or strong telluric absorption. The binned ac- 
ceptance rates of injected systems, as displayed in the 
top panel of Figure |4l motivated this functional form for 

To estimate the user acceptance rate of contaminating 
false positives, we chose a triangle function for P FP , 



P FP ( S ) EE 



pFP 

max 



p 



FP 



S < S r , 



s > s E 



(2) 



where P FP X (MLE estimate: .0474) is the maximum con- 
tamination rate, which occurs at a SNR of s p = 4.35. By 
SNR Sf ee 10, the user-error returns to zero again. The 
property that the false positive acceptance rate approach 
zero at both s = (the user never accepts a system 
with low SNR) and large signal-to-noise (the user rarely 



makes mistakes at high SNR) in conjunction with the 
lack of good motivation for a more complicated form led 
to the functional form in Equation [2j The bottom panel 
in Figure 2] displays these binned acceptance rates and 
the overplotted maximum likelihood fits. Since P Ms 11 
and Pfp are estimates of the success rate p of a binomial 
distribution, we use the Wilson Score interval to deter- 
mine the uncertainty in p, shown as error bars in each 
bin. The large errors in the bottom panel at high SNR re- 
flect a paucity of high-SNR false positives in the sample. 
While the formal uncertainty in p is large, such strong 
false positives are rare and were correctly identified in 
the few cases where they arose. 

For each individual sightline, the total acceptance rates 
pMg ii anc j pFP were estimated as a function of W r using 
the continuum-normalized error arrays to calculate the 
SNR. We denote the resulting grids of user-acceptance 
rates as A M ^ n (z,W) and A FP (z,W). The total com- 
pleteness C q (z,W) for each QSO q is then the product 
of the fraction of systems identified by the line finder 
Lq(z,W) and the fraction of systems passing visual in- 
spection, 



C q (z,W) = L q (z,W)Af^(z,W). 



3.2.3. Survey Completeness 



(3) 



The total pathlength-weighted completeness across all 
sightlines is shown in Figure [5j Figure [B] shows complete- 
ness calculations at three fiducial values of W across the 
full survey and also for the three QSO spectra depicted in 
Figure [1] The acceptance rates for false positives do not 
enter our simulated completeness results here, but will be 
combined with the rejected candidate distribution later 
in Section 231 in order to adjust the population distribu- 
tion directly. 

We can use our completeness calculations for each 
QSO C q (z,W) to calculate the redshift path density 
g(z, W) of our survey, defined as the total number of 
sightlines at a given redshift z for which Mg n sys- 
tems with rest equivalent widths greater than W are 
observable. Typically, this quantity is defined with a 
hard cutoff (|Steidel fc SargentlH992TlNestor et~all 120051: 
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Figure 6. Sightline-specific completeness for three fiducial equiv- 
alent widths W r '. 0.3, 0.5, and l.oA (from bottom to top in each 
panel). The first three panels correspond to the three QSO spectra 
shown in Figure [l] and typify the best, median, and a poor SNR 
in our sample, respectively. The fourth panel is the average across 
the full survey. 
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Redshift 

Figure 7. The completeness-weighted number of sightlines over 
which we are sensitive at three specific rest equivalent widths: 
0.3,0.5, and l.oA (from bottom to top in each panel). The top 
panel corresponds to the 35 quasars in our survey with redshift less 
than 5, the middle panel to the 11 quasars with redshift greater 
than 5, and the bottom panel to the full survey. 



iProchter et al.ll2006[ ); mathematically, 
Q 

g{z, W) = R q (z)9(W - W™ m {z)), hard cutoff (4) 

9=1 

where the sum is over QSOs, Q is the total number of 
sightlines probed, 9{x) is the Heaviside function, and 

W^ am (z) is the rest equivalent width threshold floor. 
The factor R q {z) represents the Mg II search region 
probed for QSO q, and is 1 for every value of z for 
which a Mg n system of any size could conceivably have 
been found, and everywhere else. (More specifically, 
it is 1 at every redshift within the minimum and maxi- 
mum redshift Mg n search limits for QSO q, except for 
those regions excluded because of poor telluric correc- 
tions.) The threshold floor W^ nm (z) here sets the limit 
on discovered Mg n systems that are included in the final 
analysis: Mg n systems located at redshift z m for QSO 

q with rest equivalent widths W m > W^ clm (z m ) are in- 
cluded, while smaller absorption systems are not. These 
threshold floors are specifically chosen such that the QSO 
spectra have (at least roughly) 100% completeness above 
these levels. 

If completeness is explicitly quantified, however, then it 
may be folded directly into the pathlength calculations. 
For example, having two QSO sightlines that are 50% 
complete at a given redshift is as good as having one 
which is 100% complete. If known, this completeness 
may be folded into our expression for the redshift path 
density, 

Q 

g(z,W) = J2R q (z)C q (z,W), general (5) 

9=1 

where C q (z,W) is the user-error adjusted completeness 
from Equation [3] Comparing this equation with Equa- 
tion @] above for the redshift path density with a hard 
threshold floor reveals that such a floor amounts to ap- 
proximating the completeness as either or 1 at all red- 
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Figure 8. The total path g(W) as a function of rest frame equiv- 
alent width, as given in Equation [7] 

shifts and rest frame equivalent widths, 

C q {z,W) = {\ W>wf*{zlR q {z) = l (6) 
9V ; [0 otherwise ^ ' 

Unlike Equation 0] for g(z, W) above, this formulation in 
Equation [5] does not simply exclude regions of the data 
which, although perhaps not 100% complete, still contain 
valuable information. 

The redshift path densities g{z, W) for the QSOs in 
our survey with z < 5, z > 5, and all z are shown in the 
three panels of Figure [7] Figure [8] displays the total path 
g(W) as a function of rest frame equivalent width, found 
by integrating the redshift path density g(z, W) over the 
full redshift search range of our survey, 

g(W)= f g(z,W)dz. (7) 



3.2.4. Parameter Errors 

We also used our Monte Carlo simulation to study the 
errors in our measured line parameters. Figure |H] shows 
histograms of the differences in the injected and mea- 
sured rest frame equivalent widths (top panel) and red- 
shifts (bottom panel). The rest frame equivalent width 



10 




-0.002 



-0.001 0.000 0.001 

Redshift errors 



0.002 



Figure 9. Characterization of measurement errors from our 
Monte Carlo simulations. The top panel shows equivalent width er- 
rors calculated by comparing measured to injected values; we found 
zero mean offset and cryy = 0.055A. The bottom panel shows errors 
in measured redshift, with zero mean offset and a z = 2.74 X 10~ 4 . 
For reference, this error is comparable to the size of one FIRE pixel 
(dzP lx = 2.25 X 10 — 4 ), which is indicated with solid vertical lines. 
The dotted vertical lines represent the precision provided in Table 
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Figure 10. Histograms over rest frame equivalent width (top 
panel) and redshift (bottom panel) for the 110 Mg II systems found 
in our survey and listed in Table[2] (proximate system not included). 
These plots represent raw counts; they do not include corrections 
for incompleteness or false-positive contamination. 

plot contains the errors for over 1750 Mg n 2796A and 
2803A singlets fit by hand. The overplotted zero mean 
Gaussian has a width equal to the standard deviation 
of the sample, aw = 0.055A. The plot shows no signifi- 
cant biases. The redshift plot contains errors for the first 
200 simulated spectra for each QSO in our Monte Carlo 
completeness simulation (corresponding to over 67,500 
systems). The solid vertical lines are separated by one 

pixel (d2;P lx = 2.25e — 4), and the dotted vertical lines 
represent the precision given in Table [5] The standard 
deviation of the sample is er z = 2.74e — 4. 

4. RESULTS 
4.1. Population Statistics 

Figure [TTj] shows raw histograms of the rest frame 
equivalent width (top panel) and redshift (bottom panel) 
distribution of the sample's 110 Mg n systems, before 
corrections for incompleteness (proximate system not in- 
cluded). The rest frame equivalent widths range from 
W r = 0.08 ± O.OlA (SDSS0140-0839, z = 3.71) to 
W r = 5.58±0.07A (CFQS1509, z = 3.39). The sam- 



800 



600 



| 400 

< 

200 





1 2 3 4 5 6 
Rest frame EW (A) 

Figure 11. The velocity spreads (in km/s) plotted against the 
2796A singlet rest frame equivalent width for the 110 located sys- 
tems listed in Table [2] (proximate system not included). The 
overplotted line is the sigma-clipped best fit, Av = (125.7 ± 
6.0 km/s/A)W r (2796) + (80.2 ± 8.6 km/s). 

pie contains 5 Mg II systems with redshifts over 5, with 
a maximum value of z = 5.33 (SDSS1411, W r = 0.20A). 
Three of these systems come from SDSS1411's sightlinc, 
which contains four high redshift systems at z — 4.93, 
5.06, 5.25, and 5.33. 

Figure [TT] displays the relationship between the rest 
frame equivalent widths of the 2796A Mg n singlets 
and velocity spreads (measured to where the absorption 
troughs intersect the continuum) of this sample. These 
velocity spreads have been adjusted for FIRE's resolution 
of 50 km/s. The overplotted line is the sigma-clipped 
best fit line, given by 



Av = (125.7 ±6.0) 



W r 

IX 



(80.2 ±8.6) kms" 1 . (8) 



This slope and the strong correlation prese nt in Figure 
are consistent with previous results from lEllisonl (200 1 
at lower redshifts, 0.2 < z < 2.4. We incorporated this 
measured correlation and scatter into our Monte Carlo 
completeness simulations described in 13.2.11 to ensure 
that our injected systems resembled true systems from 
the sample. (Strong injected systems were modeled as a 
blend of narrow components with total velocity spread 
drawn from the distribution specified by Equation [8]). 

4.2. Incorporating Variable Completeness in Population 
Distributions 

When calculating the line density and equivalent width 
distributions, we first define the true number of systems 
in the universe within a binned region k centered on 
(zk, Wk) as Nk- If the total redshift path in bin k is Azk, 
and the total range of rest frame equivalent width in the 
bin is AWk, then the population densities are given by 

d 2 N \ N k 
dzdW) k Az k AW k 
'dN\ = N k 
dz J k Az k 
fdN_\ JV^_ 
\dWj. AW k 



(9) 
(10) 

(11) 



Using our completeness matrices we may relate these 
quantities to the expected number of survey detections. 
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A complete derivation of this method is given in the Ap- 
pendix. In summary, for an incomplete survey with path 
length given by g(z, W), the number of detected systems 
in a bin k is estimated as 



N k (z,W) 



k 9{z > w) ^dW 



dzdW. 



(12) 



The average completeness C k is defined as the ratio be- 
tween this and the true number of systems N k , 



C k = N k /N k 



ff k {EUR q (z)c q (z,w) 



d - dzdW 

dzdW u/ -' avv 



lit (Y£,i *,(*)) £MvdzdW 



(13) 



(14) 



Our Monte Carlo calculations give us access to C k 
(which in turn leads to N k ) with an important caveat. 
The average completeness within a bin is a weighted inte- 
gral of the fine grained completeness across the bin, with 
weighting function d 2 N / dzdW . If the true completeness 
varies within a bin, then that average depends on the 
very function we are trying to calculate. 

There are three ways around this dilemma: (1) We 
can assume that the completeness C q (z,W) is constant 
across all sightlines and across the full breadth of the bin. 
Since the quality of the data varies largely from spectrum 
to spectrum, this assumption is not warranted. (2) We 
can assume that the population density within the bin is 
approximately constant. With this approximation, the 
population density in bin k reduces to 



d 2 N \ 



dzdW 



) k J k g(z,W)dzdW 



(15) 



T his is the met h od em ployed in the recent C IV study 
of iSimcoe et al.l ()201lD when considering variable com- 
pleteness within bins, and implicitly in any previous 
Mg II studies which employed the hard threshold com- 
pleteness floor given previously in Equation [6] (3) We 
can compensate for this variability in both completeness 
and absorber line density by employing a maximum like- 
lihood estimate to a functional form of the frequency dis- 
tribution in the calculation of the average completeness 
in Equation 1141 

We explored methods (2) and (3) using the maximum- 
likelihood estimates of the frequency distribution defined 
later in this paper. Figure Q21 shows the effect of these 
corrections on dN/dz for the range W > 0.3 A. "X"- 
shaped points indicate our dN/dz values prior to com- 
pleteness correction; the triangles show completeness- 
corrected values with no user-screening (i.e., using L 
rather than C for the correction and not adjusting for 
false positives, as described next section). The squares 
and diamonds include the full completeness and false pos- 
itive corrections using methods (2) and (3) respectively 
to treat intra-bin variation in the frequency distribution. 
Apparently this effect contributes a negligible correction 
to our calculations of the population densities when com- 
pared with the sample's Poisson errors. For this reason, 
we will instead use method (2) for calculating C unless 
otherwise stated. 



2 
is 



W, > 0.3 



Redshift 

Figure 12. Illustration of the effect of incompleteness corrections 
for systems with W > 0.3A. "X"-shaped points show the line den- 
sity prior to completeness correction. Triangles incorporate the au- 
tomated completeness L q (z,W) but no adjustment for interactive 
rejection of candidates. Squares further incorporate the additional 
adjustment for the user-r ejec tion of real systems and acceptance of 
false positives (Equation [16} . Diamonds incorporate initial maxi- 
mum likelihood estimates of d 2 N/dzdW to adjust for variability in 
completeness across each bin (Equation [T4} . Quoted errors reflect 
both counting statistics and uncertainty on completeness estimates, 
but the Poisson errors dominate in all cases. 

4.3. Adjusting for False Positives 

In the previous section, we calculated population den- 
sities in the presence of variable completeness, but did 
not account for false positives. Adjustments for false 
positives once again require Equations [9] to [IT] for the 
population densities, but with the true number of sys- 
tems N k calculated as 



N k = 



N k {l 



■4) 



- 4h 



Cfc 



LkA F k 



(16) 



where N k is the number of detected Mg II systems in 
bin k, F k is the number of rejected candidates, C k is the 
average completeness, L k is the automated line identifi- 

F 

cation finding probability, and A k is the user acceptance 
rate of non-Mg II candidates. All average values here are 
path length weighted (once again employing the assump- 
tion that the population distributions are approximately 
constant across the bin). A full derivation of this formula 

p 

is given in the Appendix. Note that as A k — > (the user 
correctly rejects all false positives), this reduces to the 

previous formula, C k — N k /N k . 

Figure [12] shows the effect of user errors on the resul- 
tant dN/dz for the rest frame equivalent width range 
W > 0.3A. As previously stated, the triangles employ 
neither Equation [16] above to adjust for accepted false 
positives, nor the full completeness C q (z,W), which ac- 
counts for user rejection of real Mg 11 systems. The 
squares include both of these user error corrections. The 
rejection of real Mg 11 systems is by far the more domi- 
nant of these two effects, with the correction for accept- 
ing false positives negligible compared to the Poisson 
error bars (indicative of conservative users hesitant to 
accept anything suspect). The correction is the largest 
in the lowest redshift bin, where the effective sensitivity 
is 0.393 A (as shown later), the poorest regime in this 
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1.00 



0.10 r 



0.01 r 



2 3 4 

Rest frame EW (A) 

Figure 13. The Mg II equivalent width distribution d 2 N/dzdW 
for the full survey, 1.9 < z < 6.3. The overplotted line represents 
the maximum likelihood fit to the exponential form of Equation 
1171 with best-fit parameters given in Table [3] 



Table 3 

Mg II Equivalent Width Distribution, 
Full Sample and Redshift Cuts 



(W r ) 


AW r 


C 


Number 


d 2 N/dzdW 


(A) 


(A) 


(%) 








2 


= 1.90 


- 6.30 




0.42 


0.05-0.64 


50.0 


47 


1.570±0.272 


0.94 


0.64-1.23 


81.2 


28 


0.594±0.119 


1.52 


1.23-1.82 


82.6 


14 


0.291±0.080 


2.11 


1.82-2.41 


82.6 


9 


0.187±0.064 


2.70 


2.41-3.00 


82.6 


1 


0.083±0.042 


4.34 


3.00-5.68 


82.6 


6 


0.027±0.011 




2 


= 1.95 


- 2.98 




0.43 


0.05-0.64 


44.3 


22 


1.712±0.412 


0.94 


0.64-1.23 


77.5 


8 


0.381±0.137 


1.53 


1.23-1.82 


79.2 


7 


0.323±0.125 


2.12 


1.82-2.41 


79.2 


4 


0.187±0.094 


2.71 


2.41-3.00 


79.2 


2 


0.093±0.066 


4.34 


3.00-5.68 


79.2 


5 


0.051±0.023 




2 


= 3.15 


- 3.81 




0.41 


0.05-0.64 


61.7 


17 


1.634±0.414 


0.94 


0.64-1.23 


91.1 


11 


0.723±0.222 


1.52 


1.23-1.82 


92.0 


6 


0.391±0.161 


2.11 


1.82-2.41 


92.0 


2 


0.130±0.092 


2.70 


2.41-3.00 


92.0 


1 


0.065±0.065 


4.34 


3.00-5.68 


92.0 


1 


0.014±0.014 




2 


= 4.34 


- 5.35 




0.41 


0.05-0.64 


61.0 


7 


1.429±0.559 


0.94 


0.64-1.23 


90.7 


5 


0.694±0.313 


1.52 


1.23-1.82 


91.5 


1 


0.138±0.138 


2.11 


1.82-2.41 


91.5 


2 


0.275±0.195 


2.70 


2.41-3.00 


91.5 


1 


0.138±0.138 


4.34 


3.00-5.68 


91.5 





<0.030 



survey. As the sensitivity improves with redshift, the 
user adjustments decrease in magnitude. Across the full 
range of redshifts and equivalent widths, the user accepts 
> 90% of all real Mg n systems. The presented error bars 
include completeness errors, but at < 5%, these are sub- 
dominant to the Poisson errors. 

4.4. (PN/dzdW 

Figure [T3] displays the rest frame equivalent width dis- 
tribution d 2 N/dzdW for the entire survey with all ap- 
propriate corrections applied. Table |3] contains these cal- 
culated values. The distribution is fit well by a simple 
exponential at W < 3.0A, but shows a relative overabun- 
dance of systems in the highest equivalent width bin. 



Table 4 

Maximum-Likelihood Fit Parameters for 
Exponential Parameterization of the W r 
Distribution 





A2 


W* 


N* 






(A) 





0.68 a 
1.10 a 
1.60 a 
2.51 
3.46 
4.78 
3.24 



0.366-0.871 
0.871-1.311 
1.311-2.269 
1.947-2.975 
3.150-3.805 
4.345-5.350 
1.947-6.207 



0.585±0.024 
0.741±0.032 
0.804±0.034 
0.935±0.150 
0.766±0.152 
0.700±0.180 
0.824±0.090 



1.216±0.124 
1.171±0.083 
1.267±0.092 
1.755±0.081 
1.952±0.105 
1.811±0.141 
1.827±0.059 



a Parameter fits from Nes tor et all \200S ) 

Using maximum-likelihood techniques, we fit the 
equivalent width distribution using the functional form 

d 2 N 



dzdW 



-w/w, 



(17) 



We first calculated a maximum likelihood estimate of W* , 
and then calculated A* by insisting that the pathlength- 
weighted integral of this population distribution equal 
the number of found systems, 



g(z,W) 



d 2 N 
dzdW 



-dzdW. 



(18) 



This process yielded values of W* = 0.824 ± 0.090A 
and N, = 1.827 ± 0.059. Figure HI displays this maxi- 
mum likelihood fit plotted over the binned, completeness- 
corrected values. 

We can use this result to calculate an effective equiv- 
alent width sensitivity for the full survey. Suppose one 
assumes that the survey's limiting equivalent width is 
characterized by a single number that is constant at all 
redshifts and along all sightlines, Wfl oor - Then according 
to Equation 01 

Q 

g(z,W)=6(W -VF floor )^i? g ( 2 ). (19) 

9=1 

In this regime, the maximum likelihood estimate of W* 
may be solved analytically as 



\ m=l / 



Wflo 



■w- w flc 



(20) 
(21) 



where W is the average rest fra me equivalent width of our 
sample (jMurdoch et alJll986D . We define the effective 
sensitivity W c s of our survey as the value of Wfloor in this 
equation that, given our calculated value of W*, would 
lead to the average rest frame equivalent width W that 
we observe. For our reported values of and W, we 
obtain an effective sensitivity of W e ft = 0.337 ± 0.090A 
across the full survey. 

In Figure [I4j we subdivide d 2 N/dzdW into three dif- 
ferent redshift ranges; 1.95 < z < 2.98, 3.15 < z < 3.81, 
and 4.34 < z < 5.35. Table|3]lists the values plotted. We 
also refit MLE estimates of W* for each of these regions 
(0.935 ± 0.150, 0.766 ± 0.152, and 0.700 ± 0.180, respec- 
tively), and recalculated the effective rest frame equiv- 
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0.01 
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1.95 < z < 2.98 




3.15 <z< 3.81 



4.34 < z < 5.35 







1 2 3 4 5 
Rest frame EW (A) 



Figure 14. The equivalent width distribution d 2 N/dzdW broken 
into three redshift intervals. Over plotted are the MLE fits to the 
exponential form of Equation 1171 Table \3\ contains the redshift 
ranges used and best-fit parameters for each case. 




Figure 15. Evolution in steepness of the equivalent width distri- 
bution, as parametrized by its exponential scale factor W*. The 
triangles represent maximum-likelihood esti ma tes from the FIRE 
survey for the three r edshift ranges in Fig ure ll4l The crosses repre- 
sent the values from[Ncstor ct al. (2005), determined from SDSS. 
The overplotted solid line is the maximum likelihood fit for this 
survey's data from Equatio n 1221 and the dashed line is the anal- 
ogous fit from [Nestor ct al. (2005). W* evolves significantly over 
redshift, and a simple polynomial form does not adequately capture 
this evolution. 



alent width sensitivities in each region (0.393 ± 0.150, 
0.278 ± 0.152, and 0.269 ± 0.180, respectively). In Fig- 
ure [15] we p lot W*(z), an d inclu de previously published 
values from iNestor et ahl ((2005), to connect our results 
with lower redshifts. Table [4] contains a summary of all 
these values. 

Clearly the trend of increasing W* seen at lower red- 
shift does not extrapolate to earlier epochs. Rather, it 
exhibits a peak value of ~ 0.9 at z ~ 2.5 and smaller val- 
ues at both higher and lower redshifts. A larger value of 
W* corresponds to a flatter equivalent width distribution 
with a larger fractional contribution from high-Wr ab- 
sorbers. It appears that going back in time, d 2 N/dWdz 
evolves to include more strong Mg n systems until a 
maximum is reached at z ~ 2.5, after which the rela- 
tive frequency of strong systems once again declines. It 
is tempting to associate this rise and fall of the strong 
Mg II systems with the concurrent rise and fall of the 
global star formation rate; we will explore this associa- 
tion in detail below. 

To q uantify W*'s redshi ft evolution, we followed the 
lead of INestor et all (|2005D and performed a MLE fit to 



d 2 N 



dzdW w,(i + z y 



-w/w t {i+zy 



(22) 



The best fit parameters for this from were W* = 1.166 ± 
0.359A, 5 = -0.229 ± 0.220, and N* = 1.814 ± 0.057. 
The parameters W* and <5 were anti-correlated, with 
Pw,s — —0.940. In Figure [15] we overplotted this best 
fit with a sol i d line , as well as the analogous MLE fit of 
INestor et al.l (|2005l ) to his data with a dashed line. The 
results show that a simple power law in (1 + z) does not 
adequately capture W*'s redshift evolution over the full 
redshift range now available. Rather, there is 2a evi- 
dence that the overall distribution is steepening towards 
higher z. 

To examine evolution in the overall normalization of 
the frequency distribution, we also fit our data to the 
functional form 



d 2 N 



N, 



dzdW W, 



7?f(l + 2) 



P--W/W, 



(23) 



This exercise yielded W* = 0.824 ± 0.090A, /3 = 0.002 ± 
0.487, and iV* = 2.211 ± 1.533. As is usual, these multi- 
dimensional MLE fits yield large errors for small samples 
such as ours, but within these limitations we do not see 
statistically significant evidence for redshift evolution in 
the normalization of number counts. 

4.5. dN/dz and dN/dX 

Figure [Tni displays the completeness-corrected, binned 
values of dN/dz for our Mg n sample in three differ- 
ent rest frame equivalent width ranges. These ranges 
were chosen for easy comparison w ith published SPSS 
Mg II surveys at l ower redshift (jNestor et al.l 120051 : 
iProchter et al.l f2006h . which are overplotted. Table E] 
contains the values of the data points in this figure. 

We also fit the absorber line density to a simple poly- 
nomial, 

*f=N.{l + zf. (24) 
dz 

We calculated a maximum likelihood fit of /3, and then 
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Figure 16. The absorber line density dN/dz for three different 
rest frame equivalent width ranges. The circular da ta points are 
from this study, the triangle data points are from INestor et al.l 
(120051') . and th e squa re data points (bottom panel only) are from 
IProchter et a l. (2006). Table \E\ contains the values associated with 
the FIRE data. The overplotted dashed lines represe nt m aximum 
likelihood fits to the functional form given in Equation [24] over the 
range of the FIRE data only. 

calculated AT* by insisting that our redshift path density 
integrated against this absorber cross section equal the 
total number of systems located, 

9^W)-±- d -^dzdW, (25) 



AW dz 



where AW is the total rest equivalent width range 
probed. We performed this MLE fit over our full survey, 
and over the three rest frame equivalent width ranges 
from Figure [TBI Table El c ontains these MLE fits , as well 
as the analogous fits from IProchter et al.l ([2006ft . These 
fits are also overplotted with dashed lines in Figure [16] 

In addition, we calculated the comoving absorber line 
density dN/dX. These completeness-corrected, binned 
results are sho wn in Figure fl7l w ith t he correspond- 
ing re sults from INestor et al.l (|2005l ) and IProchter et afl 
(2006) once again overplotted. The dashed horizontal 
lines show the stra ight mean dNj dX of our data com- 
bined with those of INestor et al.l (120051) . As in the usual 
interpretation, a constant value of dN/dX corresponds 
to a population with no comoving evolution. 

The smaller rest frame equivalent width systems (W r < 
lA) show remarkably little evidence of cosmic evolution 
over the entire range of redshifts probed by these sur- 
veys. For the low end 0.3 < W r < 0.6 A sample the 
data are consistent (xl = 0-2) with a constant value 
of dN/dX = 0.11, and only 16% standard deviation. 
The same is true for 0.6 < W r < 1.0, where uncertainty 
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Figure 17. Evolution of the line density per comoving length 
interval dN/dX. Over plotted horizontal li nes represent the joint 
mean of the data from Nestor ct al. (2005) and this survey. The 
top two panels, representing systems weaker than W r < l.oA, are 
statistically consistent with zero evolution from z = 0.4 to z = 5.5. 
In contrast, the stronger Mg II systems (W r > lA) grow in number 
density until z ~ 3, after which the number density declines toward 
higher z. 



Table 5 



Mg II Absorption line density dN/dz 


(*) 


Az 


C Number dN/dz 
(%) 


dN/dX 


0.3A < W r < 0.6 A 


2.231 
2.727 
3.457 
4.786 


1.947-2.461 
2.461-2.975 
3.150-3.805 
4.345-5.350 


47.3 4 0.362±0.193 
69.1 8 0.486±0.176 
77.9 10 0.453±0.146 
76.9 4 0.387±0.195 


0.115±0.061 
0.141±0.Q51 
0.119±0.038 
0.088±0.045 


0.6A < Wr < 1.0A 


2.232 
2.723 
3.458 
4.779 


1.947-2.461 
2.461-2.975 
3.150-3.805 
4.345-5.350 


65.5 5 0.343±0.156 

84.6 3 0.149±0.087 
89.9 11 0.432±0.132 
89.6 5 0.415±0.187 


0.109±0.049 
0.043±0.025 
0.113±0.035 
0.095±0.043 


W r > 1.0A 


2.231 
2.722 
3.459 
4.776 


1.947-2.461 
2.461-2.975 
3.150-3.805 
4.345-5.350 


70.4 8 0.509±0.185 

87.5 12 0.578±0.170 
92.0 12 0.461±0.135 
91.5 4 0.325±0.164 


0.161±0.059 
0.168±0.050 
0.121±0.035 
0.074±0.037 



from incompleteness is even smaller. In this W r range, 
dN/dX = 0.09 with xl = 0.5 and 25% scatter around a 
constant value. In contrast, the large rest frame equiva- 
lent width systems (W r > lA) show significant evolution 
over ours and previous surveys; xt = 3-3 for a constant 
value of dN/dX over our survey and Nestor's, and in- 
creases to xl = 20 if Prochter's W r > 1.0 A data are 
further included. Table [5] contains the values of our new 
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Table 6 

Maximum-Likelihood Estimates of the 
Line Density Evolution dN/dz =N*(1 + z) fi 



(W r ) 


AW r 


Az 


13 


N* 


(A) 


(A) 








1.17 a 


1.00-1.40 


0.35-2.3 


0.99lg;» 


u - oi -o.io 


1.58 a 


1.40-1.80 


0.35-2.3 


l-56±°;|f 


0.020i°;°5 


1.63 a 


1.00+ 


0.35-2.3 


1 40+ ' 16 


o.ost^ 


2.08 a 


1.40+ 


0.35-2.3 


1 74+O.22 
i -' 4 -0.22 


0.036l°;°« 


2.52 a 


1.80+ 


0.35-2.3 


i qn+0.30 
1 -y^_0.32 


0.0l6lg;°l 


0.824 


Full Range 


1.9-6.3 


0.133 + 0.477 


1.113 + 0.759 


0.45 


0.30-0.60 


1.9-6.3 


-0.104+0.937 


0.482+0.644 


0.79 


0.60-1.00 


1.9-6.3 


0.928+0.889 


0.093+0.121 


1.82 


1.00+ 


1.9-6.3 


-0.746+0.857 


1.301+1.555 



a Parameter fits from Prochtcr ct al. (2006), with corresponding 
upper and lower 95% confidence intervals. This survey's results 
include la errors. 

data points from this figure. 



5. DISCUSSION 

The initial results of our 46 sightline survey demon- 
strate that Mg II absorbers were already commonplace 
at z > 5. Rather remarkably, the comoving number den- 
sity of Mg II absorbers with W r < 1 A is completely fiat 
at dN/dX ^0.1 with scatter of only ~ 25% between 
z = 0.4 and z — 5, a period of > 8 Gyr. In contrast, 
the stronger Mg II absorbers grow in number density by 
a factor of ~ 2 — 3 from the present towards z ~ 3, af- 
ter which they drop off by a comparable amount toward 
z ~ 5. The differential nature of the evolution between 
strong and weak systems is reflected in the equivalent 
width distribution, which is more shallow (higher W*) at 
z ~ 3 than it is at lower and higher redshifts. The change 
in total number counts is still small since the strong sys- 
tems are subdominant at all redshifts. 

Local surveys of absorber-galaxy pairs typically find 
Mg II systems in the extended ~ 100 kpc haloes of galax- 
ies with L ~ L* or slightly below (Bergeron & Boissc 
1991 1 : iSteidel et al.lll994UChen et alJl2010aHSteidel et al.l 
200a Kacprzak et al.ll2011l ). According to the CDM pic- 
ture of structure formation, galaxies resembling present- 
day L* systems should be rare at z ~ 4 — 5. Assuming 
that Mg II absorbers are still found near galaxies at these 
epochs, it follows that the typical galaxy giving rise to 
Mg II absorption in the early universe must look quite 
different from the systems studied locally. 

The constancy of dN/dX for weak absorbers implies 
that the product na of the comoving number density and 
physical cross section of absorbing regions also remains 
constant. A simple interpretation is that these absorbers 
arise near the assembling progenitors of present-day L* 
galaxies. This would require that galaxy haloes were 
populated with Mg II very early in their history, and 
that the absorption properties were established before 
their stellar and ISM components were fully formed. 

A detailed analysis of absorption in individual 
Mg II systems will be presented in a companion paper; 
however in a qualitative sense it appears that Mg n sys- 
tems at z ~ 5 are very similar to those at low redshift, in 
contrast to the galaxy population. This either requires 
the absorption systems to persist in a steady state over 
much of the Hubble time, or else be replenished peri- 
odically. If Mg II systems are replenished through stel- 



lar feedback, this process must proceed in a way that 
does not produce a cumulative growth of the gas radius 
and/or filling factor of Mg n . This would seem to be a 
challenge unless some porti on of the gas re-accretes a nd 
cycles back into the galaxy (jOppenheimer et aT1l2010|) or 
else the Mg n -bearing clouds are out of equilibrium and 
transition to lower densities and higher ionization states. 

5.1. Mg II and Cold Accretion Flows 

Recently Mg II has been suggested as a candidate 
tracer of "cold" accretion onto galaxies, since its rel- 
atively low ionization potential (1.1 Ryd) does not al- 
low for large Mg II fractions at high T. Cold accretion 
is thought to be the dominant mode of gas transport 
onto galactic disks; it is particularly e fficient in the early 
universe for sma ll , growing haloe s (Kercs et all 120051: 
IDekel et al] 120091: iFaucher-Giguere fc Keresl I2010D and 
could provide another "renewable" source for cool ab- 
sorbers. While attractive, this model for the Mg n pop- 
ulation faces some complications at the high redshifts 
probed by our survey. Principally, at these epochs, even 
though cold accre tion is much stronger (M oc (1 + z) 2 25 ; 
IDekel et al] 120091) and may present a larger absorption 
cross section, the heavy element content of the IGM 
gas reservoir is quite low: Z/Zq ~ 3 x 10~ 4 at z ~ 4 
(|Simcoel 120111 ) and even lower at earlier times. In gen- 
eral, cold stre ams may manifest as metal-poor Lyman 
limit systems (|Fumagalli et al.l 1201 11) more prominently 
than Mg II systems. We may already be seeing this phe- 
nomenon, in that the density of Lyman limit systems and 
Mg II are of comparable magnitude at low redshift, but 
unlike M g II , the LLS density continues increasing with 
redshift (jProchaska et aJJl201Clt IRibaudo et aLll2011l ) be- 
y ond z > 4. 

iTinker fc Chenl (|2010D present a model of the evolving 
Mg II population which maps Mg n absorbers onto the 
dark matter halo mass function, as a way of studying evo- 
lution of the underlying galaxies. In their model, most 
Mg II absorbers arise in haloes with M ~ 10 11 - 10 12 M© 
which are presumed to be growing via cold accretion. 
Shock heating reduces the supply of Mg n in larger mass 
haloes, an effect which is required to match the observed 
anti-correlatio n between galaxy mass and W r in the lo - 
cal universe (jBouche et al.l l200l iGauthier et al] 120091) . 
However, an additional suppression of absorption effi- 
ciency is also needed in low mass (< 10 11 M©) haloes to 
offset the rapidly rising halo mass function in this regime. 

This low-mass cutoff halts the overproduction of weak 
Mg II systems in low mass haloes, aligning the model 
with observations at z ~ 0.5. However at high red- 
shift, a larger portion of the mass is concentrated in 
these very haloes, most of which fall below the low-mass 
cut for harboring Mg n absorption. As a direct conse- 
quence, the simple halo-occupation model under-predicts 
the measured dN/dX at z > 2.5 by orders of magnitude, 
as seen in Figure [18] It may be possible to reconcile 
the model by boosting the absorption efficiency of small 
mass haloes at early times, or else growing the absorp- 
tion cross-section per halo from 1/3 of the virial radius 
at z — 0.5 to > 2R V i r at early times. However these so- 
lutions are not physically motivated a priori and would 
need further exploration. 

It appears Mg II will be a difficult ion to use as a tracer 
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Figure 18. The linear density dN/dX for strong Mg II systems 
<W r > 1.0A) from IProchter et al.l pOM) (squares'). [Nestor et aO 
(2005) (triangles) and this survey (circles), compared with the lin- 
ear density dN/d X derived from t he ha lo occupation distribution 
(HOD) model of ITinker k Ohenl 12010) (solid line). The HOD 
model under-predicts the observed density at z > 3 because of 
a paucity of large haloes at early times. To narrow this discrep- 
ancy, the absorption efficiency of low-mass haloes would need to 
be enhanced in the early universe. 



Rq. According to the studies cited above, L* g evolves 
with redshift toward brighter values in the past. With 
a constant value of Lq defining the halo scale, this im- 
plies that galaxies at L* B will have gas haloes increasingly 
larger than Rq (i.e. > 75ft. -1 kpc) in the more distant 
past. For this reason, we will find it useful to define 
the term R*(z), which denotes the gas radius of an L* B 
galaxy at any redshift. 

In the second approach, we associate gas haloes of 
size Rq with galaxies of L = L* B (z) from the luminos- 
ity function at each epoch. This is equivalent to setting 
R*(z) = R and replacing Lq with L* B (z) in the notation 
above. Since L* B increases with z, in this scenario galax- 
ies at a fixed luminosity would have smaller gas haloes 
in the past than in the present day. 

We next model the absorbers' covering fraction n as 
a function of galaxy impact parameter p, u sing an em- 
pirical approximation based on the results of iChen et al.l 
(|2010aft : 



k(p, L Bl z) 



P-R(Lb) 



a(L t 





-R(Le 



p < a(L B ) 

a<p<R(L B ) (27) 
p>R(L B ). 



of cold accretion flows. The competition between declin- 
ing IGM metal abundance and rising cold accretion cross 
section at earlier times would require a fine balance to 
achieve the very flat trend we observe in dN/dX. While 
this solution is possible in principle, it will be difficult to 
interpret. 

5.2. Comparison with Known Galaxy Populations 

When comparing our Mg n number counts with galaxy 
populations, another approach is to skip the halo mass 
function altogether, and apply simple empirical scalings 
between Mg II gas halo size, covering fraction, and galaxy 
luminosity to predict dN/dX. The model cross sections 
are integrated over measured luminosity functions as a 
function of redshift. 

To es timate each galaxy' s cross section, we use the scal- 
ings of IChen et al.l (pOlOaD . who studied a spectroscopic 
sample of 94 galaxy-absorber pairs at (z) — 0.2. Chen 
found that galaxies with higher B band luminosities pos- 
sess more extende d Mg n absorbing haloes (see also 
iSteidel et al.lll~9 94), and fit the extent of these gaseous 
haloes R to a Holmberg-like scaling, 







(26) 



where Rq 



75ft" 1 kpc, /3 = 0.35 ± 0.03, and L* 



is the luminosity associated with objects of M B = 
— 19.8 4- 51og(ft) at their survey redshift. We use Equa- 
tion [21] together with measurements of t he B band lumi- 
nosity function a t increasing red s hift ( Poli et ah! 120031: 
Wolf et all [200l iGabasch et al.l 120041 iWillmer et al l 
2006HMarchesini et al.ll2007| ) to calculate the abundance 
and cross section of absorbers and compare with our mea- 
surements of dN/dX. 

When extrapolating Equation [51] to higher redshift, we 
explored two ways of evolving the scaling. In the simplest 
approach, Rq and Lq are fixed at all redshifts, so galaxies 
with M B = —19.8 + 5 log ft always have haloes of radius 



Briefly, each galaxy has a central core of radius a(L B ) 
within which k — 100%; between a (the core radius) 
and R (the halo boundary) the covering fraction declines 
linearly to zero. We assume a mimics i?'s scaling with 
L B hi both scenarios considered, but with ao = 30ft -1 
kpc. 

The total absorption cross section a(L B ,z) presented 
by a galaxy is then 

rR(L B ,z) 

a(L B ,z)= 2npK(p,L B ,z)dp (28) 

Jo 



L* B (z) 



2d 



a*(z). 



(29) 



Here a* (z) is the cross section of an L* B galaxy's gas halo, 
which may evolve with redshift but does not depend on 
luminosity, 



o-*{z) = R,itR*(z) 2 



3 \R Rq 



(30) 
(31) 



The quantity R is the average covering fraction of the 
gaseous halo, and is R = 0.52 with the given values of ao 
and Rq. 

The Mg II frequency dN/dX may be derived by com- 
bining this per-galaxy cross section with the B band lu- 
minosity function (j){L Bl z): 

dX c 

= J a{L B ,z) <j){L B ,z) dL B , (32) 

where L m ; n is a lower luminosity cutoff, which is not 
specified a priori, but rather chosen to matc h the nor- 
maliza tion of the observed dN/dX. Using a iSchechterl 
(1976) function for 4>{L B ,z), we may calculate the linear 
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Figure 19. Line ar densities dN/dX from our survey (large, black 
circles) and from Nestor ct al. (2005) (smaller, black circles) for 
W r > 0.3A, overplotted with dN/dX derived with Equation 1331 
using rcdshift-dependent B band luminosity functions listed in the 
text. The top and bottom panels use different assumptions about 
the scaling of halo radius with luminosity: In the top panel, we 
assume the halo size Ro is a fixed value and is associated with 
galaxies of L* B at each epoch. In the bottom panel, we assume Ro 
is fixed, but associated with a fixed luminosity Lq at all redshifts. 
The smooth curves are derived from integrating redshift-dependent 
Schcchter-function fits to the luminosity function; these are solid 
over the redshift ranges where measured, and high-rcdshift extrap- 
olations are indicated with dashed lines. 



density dN/dX as 

= -^-RkR*(z) 2 ^(z)V \ n(: ) + ■> ; \- I . — — 
dX Hq \ l b\ z ) 



( 33 ) 

where T(s, x) is the upper incomplete gamma func- 
tion. This expressi on is similar to one derived in 
IChurchill et al.l (|1999f ) . the only difference being the fac- 
tor of R correction for partial covering. 

We evaluated Equati on l33l usi n g the B band lu- 
minos it y functions of iPoliet al.l 120031); iWolf et al.l 



(2003) ;fGabasch et al.l H2004D ; iWillmer et al.l |2006j); and 
iMarchesini et al.l ()2007f ) to calculate dN/dX from z = 
to z = 5. The results are shown in Figure Q1J1 with 
triangles and circles representing the predicted and mea- 
sured dN/dX (for W r > 0.3A), respectively. The solid 
curves indicate the predicted evolution based on redshift- 
depe ndent fits to th e lum inosi ty function paramete rs 
from lGiallongo et al.l (pOOlit ) and lGabasch et al.l (pOOl . 

The top panel shows results for the model where the 
luminosity associated with i?o _ sized haloes tracks L* B (z). 
In this case, the observed and predicted dN/dX match 
best when L m i n /L* B = 0.024 ±0.005. However the model 



predicts a downward evolution in dN / dX with z that is 
slightly steeper than what is observed. This means that 
the typical cross section per galaxy is higher than the 
model predicts at high redshift, suggesting that either 
the haloes are growing slightly in radius or their filling 
factors are increasing. 

The bottom panel shows the simple model where the 
fiducial luminosity associated with Rq haloes is fixed at 
all redshifts. Our best fit is obtained with L m i n /L B = 
0.2656 ± 0.0139, almost an order of magnitude larger 
than the previous model. The redshift dependence for 
this model is flatter, in better agreement with our sur- 
vey data. The relatively high value of L m i n indicates that 
a large population of very low luminosity galaxies need 
not be invoked to explain the frequency of Mg II systems 
in this model; for the most part, all W r > 0.3 Mg II sys- 
tems surround galaxies we can readily observe, even at 
high redshift. 

This picture could be invalidated if significant num- 
bers of Mg II systems are discovered near galaxies with 
L < Q.2L* B . However to date such dwarf-Mg II associ- 
ations have not been estab lished except in a few cases. 
Even in lChen et al.l (|2010aD 's sample of 94 pairs, < 10% 
of the galaxies meet this criterion. They are predomi- 
nantly at very small impact parameter (~ 10/i -1 kpc) 
and even then show only ~ 50% coverage (in contrast to 
our model, in which the coverage is 100% out to ~ 30/i -1 
kpc). 

The dN/dX points in Figure [19] apply for W r > 
0.3A since this was the range over which the R(L) corre- 
lation was first determined. However we have also seen 
that if we restrict further to W r > l.OA, dN/dX evolves 
to a peak at z ~ 3 and declines thereafter. Neither of our 
halo models is able to reproduce this feature, which indi- 
cates that other physics must be at play. If the systems 
giving rise to the strong absorbers are contained in the 
B band luminosity function, then their gas cross section 
and/or covering fraction must rise and then fall again. 
Alternatively they may reflect systems whose comoving 
number density rises and falls without being picked up 
in rest frame B band surveys. 

Rather than start with a halo and coverage model and 
derive dN/dX, we may instead start with the observed 
Mg II dN/dX and estimate the effective size of gaseous 
haloes. To do this, we define the effective halo radius 
R (z) of an absorber according to 



wR ee {z) 2 = 7r KL (z)R L (z) 2 

_H 1 dN 
c n(z, i mm ) dX 



(34) 
(35) 



where kl{z) is the luminosity- weighted covering fraction 
and Rl{z) is luminosity- weighted gaseous halo radius. 

In Figure [201 we plot the effective halo radius as a func- 
tion of redshift, normalized by the radius calculated in 
the lowest redshift bin. We use an avera ge of the redshift 
evolut ion lu minosity functi o n mod els of Giallong o et al.l 
( 2005) and iGabasch et al.l (|2004l ) in conjunction with 
four different lower luminosity cutoffs: 0.001, 0.01, 0.1, 
and 1.0 L* B (z). All four choices of a lower luminosity cut- 
off result in the effective absorbing halo radius R e ^{z) 
increasing by ~ 40% from today to z ~ 3 — 4, and then 
slightly decreasing again (although the error bars are too 
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Figure 20. Effective halo radius (Equation I35H as a function of 
redshift for t he ob served Mg II dN/dX of this study and from 
INestor ct al. (2005), normalized by its value in the lowest redshift 
bin. We derived the absorber number density n(z, £ mm ) from a 
co mbination of th e redshift depe ndent luminosity function models 
of Giallongo ct al. (2005) and Gabasch et al. (2004) for four lower 
luminosity cutoffs: 0.001, 0.01, 0.1 and 1.0 L* B (z). With all four 
choices, the effective halo radius increases with redshift by about ~ 
40% from today until z ~ 3 — 4, and then decreases again (although 
this decrease is not statistically significant given our error bars). 

large to consider this final decrease statistically signifi- 
cant). This increase in effective halo radius must result 
from either an increase in the luminosity- weighted cover- 
ing fraction kl{z) at earlier times, or an increase in the 
extent of the gaseous absorbing haloes Rl(z) at earlier 
times. Although the lower luminosity cutoff still remains 
a free parameter in this model, the strong quantitative 
agreement of the normalized effective halo radii across 
three orders of magnitude in L m ^ n suggests this feature 
is largely independent of the lower luminosity cutoff cm- 
ployed. 

5.3. Deriving the SFR Density from dN/dX 

Thus far we have focused the discussion on moderate 
strength Mg n absorbers which do not evolve with red- 
shift. But the systems with W r > lA do evolve, with a 
characteristic rise and fall that suggests a possible con- 
nection to the cosmic star formation rate density 

Thi s connection is also explored in IMenard et al.l 
(|2011D , where a method is outlined for deriving the SFR 
density from the Mg n population distribution by es- 
tablishing a connection between Mg II rest frame equiv- 
alent width W r and [O n] luminosity, which previous 
studies have shown may be used a s a gauge of the SFR 
([Gallagher et al.l 119891: iKennicuttJ Il992t iHopkins et al.l 
120031: IKewlev et all 12004 IMouhcine et al.l 120051) . Since 
the detection of Mg n using spectroscopic absorption 
does not depend on luminosity, distance, or dust extinc- 
tion, such a probe would represent a valuable asset in 
pinning down the star formation history of the universe. 

In their paper, the authors search for [O n] emission 
lines from Mg n host galaxies imprinted on a sample of 
QSO spectra containing over 8500 known Mg II absorbers 
at redshifts 0.36 < z < 1.3. They detect a 15<r correla- 
tion between the absorber strength W r and the median 
[O II] luminosity surface density {Er, OJI ), which they fit 
to the functional form 

(Sl °" ) = j4 (ix) Q ' (36) 
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Figure 21. The star formation rate density derived using the 
method of lMenard et al.l 1120111) on lar ge M g II systems (W r > lA ) 
from this survey, Nestor ct al, (2005), and Prochtcr ct al. (2006) 
(circles; larger circles are this survey). Recent star formation rate 
density calc ulations derived from obse rvations with t he Hubble 
Telescope in Bouwcns ct al. (2010) and Bouwcns ct al. (2011) are 
overplotted with triangles. The method of Menard ct al. 1(20111 ) 
produces good qualitative agreement with the Hubble results, al- 
though the error bars are quite large. Applying this Mg II conver- 
sion to weaker Mg II systems (W r < lA) does not lead to a match 
with these Hubble rates. 

where A = (1.48 ± 0.18) x 10 37 ergs s" 1 kpc 2 and 
a = 1.75 ±0.11. The authors calculate the co-moving 
[O n] luminosity density Con(z) probed by Mg n ab- 
sorbers from this by using the completeness c orrected 
population distribution from N estor et al.l (|2005h accord- 
ing to 

f d 2 N 

Coi ^ = J dW Wdx^- (37) 

Finally, they convert this luminosity density Con(z) to 
a luminosity Lou and derive a SFR density by using 
an average over scalings from lZhu et~aT1 (|2009T >. 

= 4.2e-41 x Lou M yr' 1 Mpc~ 3 . (38) 

Using thi s technique. IMenard et al.l (|2011l ) convert the 
data from lQuider et al.l fj201 lh to SFR d ensities, and the n 
overplot SFR density data provided in iHopkinsI (|2004| ) . 
The results show agreement up to z = 2. 2. 

More recently, iLopez fc Chenl (|2011l ) have argued 
that the W r -L^pir\ correlation reported in lMenard et al.l 
(|2011f ) does not reflect a causal connection between 
Mg II absorption and star formation. Rather, they 
demonstrate that the observed correlation could plau- 
sibly arise from a combination of the empirically ob- 
served decline in W r at increasing galactocentric radius, 
and differential loss of galaxy light from the finite fiber 
aperture of the SDSS spectrograph. The authors run a 
large Monte Carlo simulation which reproduces th e W r - 
(T,l oii ) relationship found in IMenard et al.l (|2011l ). 

Nevertheless, in light of the observed evolution in 
dN/dX for W r > lA, the large, likely non-gravitational 
velocity spreads observed for strong Mg n systems 
(Figure [TTj) . and the demonstrated connection be- 
tween the s trongest Mg II absorbers and star form- 
ing galaxies (iBouche et aij|2007tlNoterdaeme et al.ll2010l: 
INestor et al.ll2011l) . it is tempting to explore Menard's 
method, while keeping in mind its possible limitations. 
The broad redshift range of our sample, combined with 



Sample and Evolution of the Mg n Frequency 



19 



those of the SDSS studies, puts us in unique position to 
compare with the cosmic star formation history. 

In Figure[2l] we show the SFR density calculated using 
Equations [36] to [38] for Mg n absorbers having W r > 1 
A (large circles) . We also include prev ious results at 
low re dshift from lNestor et all ()2005D and lProchter et all 
(2006) over the same equivalent width range (smaller cir- 
cles). For comparison, we overplot (with triangles) star 
formation rate densities obtained using recent determi- 
nations based on dropout c ounts from deep HST / WFC 3 
galaxy surveys at z > 4 (jBouwens et al.l [20101 1201 1| ). 
from a combination of ground-based spect roscopic Ha 
and Spitzer MIPS 2A^m data at z ~ 3 (jReddv et al.l 
l2008h . and from the GALEX VIMOS-VLT Deep Sur- 
vey at lower redshifts (jArnouts et al.l [20051 ) . Menard's 
transformation produces star formation rate densities in 
excellent agreement with the galaxy survey literature 
throughout the entire Mg II absorption redshift search 
range probed, 0.5 < z < 5. These results do not conclu- 
sively demonstrate the accuracy of Equation [36] — indeed 
application of Menard's redshift dependent version of 
Equation [SB] led to very high SFR values per absorber at 
z ~ 5 and therefore predicted significantly higher SFR 
densities. However this exercise certainly suggests that 
the connection between Mg II absorption and SFR is real 
for the strongest subset of Mg n absorbers, and stretches 
over the full cosmic history. 

Application of the Mg II -SFR conversion to weaker 
absorbers (W r < lA) does not produce good agreement 
with galaxy surveys. These systems, as seen in FigurefTTl 
exhibit remarkably little evolution in dN/dX over all red- 
shifts considered (0.5 < z < 5). This suggests that any 
physical connection between Mg n and star formation is 
mostly limited to the stronger end of the Mg II equivalent 
width distribution. Indeed different astrophysical pro- 
cesses may govern the evolution of systems at high and 
low equivalent width, with the division occurring roughly 
at Wr ~ lA as suggested bylLovegrove fc Simcoel ()2011[ ) 
and lKacorzak fc ChurchilH (|2011h . 

6. CONCLUSIONS 

Using FIRE, we have conducted an infrared survey of 
intervening Mg II absorption toward the sightlines of 46 
distant quasars with redshifts 3.55 < z < 6.28. The over- 
all survey has an effective sensitivity of 0.337 ± 0.090A 
and a FWHM resolution of 50 km/s. A combination of 
automated and human methods were used to identify and 
measure the properties of detected systems, and we per- 
formed extensive simulations to characterize our sample's 
incompleteness and contamination from false positives. 
Our main findings are as follows: 

1. We find 110 isolated Mg II systems (and one proxi- 
mate system) ranging in rest equivalent width from 
W r = 0.08A to W r = 5.58A and in redshift from 
z = 1.98 to z = 5.33, including 5 systems with 
z > 5. 

2. The population distribution d 2 N/dzdW resembles 
an exponential with a characteristic scale W» = 
0.824 ± 0.090A across the full survey. When 
we divide the survey into three redshift regions 
and combine our results with earlier results from 



iNestor et all ([2005h . we find that this character- 
istic scale W* rises with redshift until z ~ 2 — 3 
(where we find more strong systems) , and then falls 
towards higher redshifts. 

3. For weaker Mg II absorbers (W r < 1 A), the lin- 
ear density is statistically consistent with no evo- 
lution, having dN/dX ~ 0.1 with a scatter of only 
~ 15% — 25% from z — 0.4 to z > 5, a period 
of > 8 Gyr. In contrast, the stronger Mg II ab- 
sorbers (W r > 1 A) grow in number density by a 
factor of ~ 2 — 3 from the present towards z ~ 3, 
after which they drop off by a comparable amount 
toward z ~ 5. We interpret this as evidence of 
two distinct Mg n populations, one which follows 
the global star formation rate density, and another 
which is nearly constant in comoving aggregate 
cross section. 

4. For strong Mg II absorbers (W r > 1 A), the Mg II- 
to-sta r formation rate scaling of iMenard et al.l 
fj201 lh . together with our measured Mg II fre- 
quency, agrees remarkably well with classically 
measured star for mation density r ates. While the 
recent analysis of lLopez fc Chen! (|2011l ) cautions 
against a detailed exploration of this specific model, 
the excellent agreement with our dN/dX measure- 
ments are highly suggestive of a connection be- 
tween strong Mg II absorption and star formation. 

5. For the weaker systems, models based on the pre- 
dicted Mg II occupation of dark matter haloes sub- 
stantially under-predict dN/dX at high redshifts. 
This discrepancy derives from the low absorp- 
tion efficiency of small mass haloes locally (which 
are used for scaling), combined with the relative 
scarcity of high mass haloes in the early universe. 
To resolve the difference, the numerous low-mass 
haloes present at early times would need to have an 
enhanced absorption efficiency relative to present 
day levels. 

6. We are able to produce the observed dN/dX more 
accurately using scaling relations between galaxy 
luminosity and halo size derived at low redshift, 
but integrated over appropriate high redshift de- 
terminations of the galaxy luminosity function. 
The correspondence is most accurate when the 
Mg II absorbing halo's size scales with a fixed, 
redshift-independent fiducial luminosity. Since L* B 
increases with redshift, this implies that the effec- 
tive gaseous halo of an L* B galaxy is larger at earlier 
times. This model leaves the lower integration limit 
Lmin down the faint-end slope as a free parameter, 
set to match the normalization of dN/dX. Our 
best match is for L m in ~ 0.2L|j, so that large pop- 
ulations of low-luminosity galaxies are not strictly 
required to explain the observed Mg II frequency, 
even at high redshift. 

Our results indicate that Mg n absorption remains 
commonplace in the early universe, having evolved rel- 
atively little over a period exceeding 8 Gyr. This is 
in contrast to the C iv number counts, which decline 
by almost an order of magnitude over the same interval 
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(K. Cooksey, private communication), and the density of 
the intergalactic H I Lyman alpha forest, w hich thickens 
mark edly toward the epoch of reionization ( Beck er~et al.l 
l2007h . 

At low redshift, Mg II absorbers are typically associ- 
ated with ~ 100 kpc gas haloes surrounding ~ L* galax- 
ies. The Mg II counts at high redshift can also plausibly 
be reproduced using fairly luminous galaxy populations. 
As will be shown in forthcoming work, the H I and other 
metal line properties of the high-redshift Mg II systems 
also do not differ dramatically from the local population. 

Yet the properties of the host galaxies themselves must 
have evolved substantially over th e same period. As 
estimated by iGabasch et al.l (|2004D L* B brightens by a 
factor of ~ 2, while declines by an order of mag- 
nitude over the redshift path where Mg II has now 
been characterized. If the high redshift Mg II ab- 
sorbers are similarly associated with distinct galaxies, 
then the mass substructure within their dark matter 
haloes must have evolved significantly from z ~ 5 to 
the present. Moreover, high resolution, rest-frame opti- 
cal morphology studies of distant galaxies generally find 
that at fixed stellar mass, galaxies have smaller half light 
radii (0.7 — 3 kpc at z ~ 3) going b ackwards in time 
llLaw et alJl201lUPapovich et al.l l2005; Fr anx et alj|200l 
Ivan Dokkum et al.ll2010D . 

It appears that as the stellar populations of typical 
galaxies grow from the inside out, the gaseous haloes of 
the same galaxies may be populated with metals very 
early on, and that the gas halo properties evolve much 
more weakly than do the galaxies themselves. If our 
hypothesis is correct that the observed dN/dz may be 
recovered entirely by galaxies with L > 0.2L*, then it 
should be possible to observe such systems directly at 
intermediate redshifts, up to z ~ 3 — 4. Along with 
further spectroscopy to improve upon the modest sample 
sizes presented here, such follow up of Mg II -selected 
galaxies may constitute a fruitful means for studying the 
concurrent buildup of the stellar and near field circum- 
galactic environments. 
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APPENDIX 

MATHEMATICAL FRAMEWORK FOR VARIABLE 
COMPLETENESS 

If our data set were 100% complete over the Mg II red- 
shift search range and rest frame equivalent width range 

W™ m < W < spanned by a bin A:, then the ex- 

pected number of Mg n systems Nk found within this 
bin would be given by the integral of the population dis- 
tribution d 2 N/dzdW weighted by the sum of the regions 



R q {z) probed along each sightline; more specifically, 



N k = 




\ d 2 N 

sznndzw, cm. 

q=l 



dzdW 



where the range of the integral, denoted by subscript k, 
is over the redshift and rest frame equivalent width range 
of the bin. The total path Azk in this bin would then be 



Az k 




q=l 



R q {z) dz, 



(A2) 



where the integral is now only over redshift. If we denote 
the change in equivalent width of the bin as AWk = 



as 



max 



Wfc j then the population density is estimated 



/ d 2 N 



[dzdWJk Az k AW k 



(A3) 



With less than 100% completeness, the expected num- 
ber of Mg II systems Nk found within this bin is instead 



N k = 



J l(±R q ( Z )C q (z,W^ ^dzdWm 



k^^V 



dzdW. 



(A5) 



We define the proportionality constant between Nk and 
Nk to be the average completeness Ck in this bin, 

N k = C k N k , (A6) 

which leads to 

_ _ II k (T,UR q (z)c q (z,wj) ^dzdW 



c k 



(A7) 



Mathematically, this is equivalent to calculating the av- 
erage completeness across the bin using a probability 
density proportional to the product of the population 
distribution and the summed regions probed. 

According to Equations IA3I and IA6I above, the 
completeness-corrected estimate of the population den- 
sity is then 



/ d 2 N \ N k 
\dzdW) k C k Az k AWk 



(A8) 



Similarly, the estimates of the rest frame equivalent 
width distribution dN/dW and the absorber line density 
dN / dz are given by 



dN 



N k 



and 



dN 
~dZ 



Nk 



C k Az k 



(A9) 



(A10) 



If, as we argue in the main text, we assume that the 
population distribution d 2 N jdzdW is constant over bin 
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k when calculating Ck, then the average completeness 
across the bin reduces to 



_ / L (Eli Rq(z)C q (z, W) ) dzdW 



of a single bin. Therefore, we may define the expected 
value of an arbitrary function h(z, W) over bin k accord- 
ing to 



Sk{T.%iR^))dzdw 

J J k g(z,W)dzdW 
Az k AW k 



(All) 



(A12) 



where we have employed the definition of g(z, W) from 
Equation [5] and the definition of the total path length 
given in Equation IA2I Plugging this equation for Ck 
back into Equation IA8I for the population density above 
yields 

f d 2 N \ = N k 



\dzdW J k J k g(z,W)dzdW 



(A13) 



MATHEMATICAL FRAMEWORK FOR FALSE 
POSITIVE CORRECTION 

Let d 2 F I dzdW denote the population density of false 
positives that our automated line identification system 
flags as potential candidates. The total number of false 
positives Fk in bin k is given by 




R q (z) 



9=1 



d 2 F 
dzdW 



-dzdW. 



(A14) 



If Ay (z, W) is the probability of the user accepting a 
false positive in sightline q at redshift z and rest frame 
equivalent width W, then the number of actual false pos- 
itives F k that the user correctly identifies is 




Q 

5> 9 (z) (l-A*(z,W)] 

9=1 



d 2 F 
dzdW 



dzdW. 



(A15) 

Along with these, the user incorrectly identifies some real 
Mg ii systems as false positives, 



d 2 N 
dzdW 



P k= J J^fe^W {L q {z,W)-C q {z,W))^ 

(A16) 

where L q (z,W) is the automated completeness and 
C q (z, W) is the total completeness as given in Equation 
02 which also accounts for the user's decisions. 

In terms of accepted candidates, the number of cor- 
rect ly ac cepted Mg II candidates N k follows from Equa- 
tion El 



« \ d 2 N 

R q (z)C q (z,W)\ j^dzdW, (A17) 



dzdW, 



Nt = 




Along with these systems, the user incorrectly accepts 
some false positives, the number N k of which is given by 




Q n 
Y,R q {z)A F q {z,W) 

9=1 / 



d 2 F 
dzdW 



dzdW. (A18) 



As before, we will assume that the Mg n and false pos- 
itive population densities vary negligibly over the course 



hk = 



I k\T. Q q =iR q {z)h q {z J W)) dzdW 



(A19) 



With this definition, the number of identified (real and 
spurious) Mg II systems Nk and false positives Fk in bin 
k are given by 

N k = N° + Nl = NkCk + FkA^ 

F k = F° + Fl = (1 - Af)Fk + N k (L k - C k ). 

Here, we have two equations and two unknowns, namely 
the number of real systems Nk and true false positives 
Fk- Solving for Nk gives 



Nkjl-A^-A^Fk 

— F 

C\ — LkA k 



(A20) 



The population densities are then given by Equations [9] 
to [TT] with this value substituted for Nk- 

REFERENCES 
Arnouts, S., et al. 2005, ApJ, 619, L43 

Becker, G. D., Rauch, M., & Sargent, W. L. W. 2007, ApJ, 662, 
72 

Bergeron, J., & Boisse, P. 1991, A&A, 243, 344 

Bochanski, J. J., et al. 2009, PASP, 121, 1409 

Bouche, N., Murphy, M. T., Peroux, C, Csabai, I., & Wild, V. 

2006, MNRAS, 371, 495 
Bouche, N., Murphy, M. T., Peroux, C, Davies, R., Eisenhauer, 

F., Forster Schreiber, N. M., & Tacconi, L. 2007, ApJ, 669, L5 
Bouwens, R. J., et al. 2010, ApJ, 725, 1587 
— . 2011, ApJ, 737, 90 

Bowen, D. V., k, Chelouche, D. 2011, ApJ, 727, 47 

Chen, H.-W., Helsby, J. E., Gauthier, J.-R., Shectman, S. A., 

Thompson, I. B., & Tinker, J. L. 2010a, ApJ, 714, 1521 
Chen, H.-W., Wild, V., Tinker, J. L., Gauthier, J.-R., Helsby, 

J. E., Shectman, S. A., & Thompson, I. B. 2010b, ApJ, 724, 

L176 

Churchill, C. W., Rigby, J. R., Charlton, J. C, & Vogt, S. S. 

1999, ApJS, 120, 51 
Cushing, M. C, Vacca, W. D., & Rayner, J. T. 2004, PASP, 116, 

362 

Dekel, A., et al. 2009, Nature, 457, 451 
Ellison, S. L. 2006, MNRAS, 368, 335 

Elston, R., Bechtold, J., Hill, G. J., & Ge, J. 1996, ApJ, 456, L13 

Faucher-Giguere, C, & Keres, D. 2010, ArXiv e-prints 

Franx, M., van Dokkum, P. G, Schreiber, N. M. F., Wuyts, S., 

Labbe, I., & Toft, S. 2008, ApJ, 688, 770 
Fumagalli, M., O'Meara, J. M., & Prochaska, J. X. 2011, Science, 

334, 1245 

Gabasch, A., et al. 2004, A&A, 421, 41 

Gallagher, J. S., Hunter, D. A., & Bushouse, H. 1989, AJ, 97, 700 
Gauthier, J., Chen, H., & Tinker, J. L. 2009, ApJ, 702, 50 
Gauthier, J.-R., & Chen, H.-W. 2011, MNRAS, 1604 
Gauthier, J.-R., Chen, H.-W., & Tinker, J. L. 2010, ApJ, 716, 
1263 

Giallongo, E., Salimbeni, S., Menci, N., Zamorani, G., Fontana, 
A., Dickinson, M., Cristiani, S., & Pozzetti, L. 2005, ApJ, 622, 
116 

Hopkins, A. M. 2004, ApJ, 615, 209 

Hopkins, A. M., et al. 2003, ApJ, 599, 971 

Jiang, L., Fan, X., Vestergaard, M., Kurk, J. D., Walter, F., 

Kelly, B. C, & Strauss, M. A. 2007, AJ, 134, 1150 
Kacprzak, G. G, & Churchill, C. W. 2011, ArXiv e-prints 



22 



Kacprzak, G. G., Churchill, C. W., Evans, J. L., Murphy, M. T., 

& Steidel, C. C. 2011, MNRAS, 416, 3118 
Kennicutt, Jr., R. C. 1992, ApJ, 388, 310 

Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, MNRAS, 
363, 2 

Kewley, L. J., Geller, M. J., & Jansen, R. A. 2004, AJ, 127, 2002 
Kondo, S., ct al. 2008, in Astronomical Society of the Pacific 
Conference Series, Vol. 399, Panoramic Views of Galaxy 
Formation and Evolution, ed. T. Kodama, T. Yamada, & 
K. Aoki, 209 

Lanzetta, K. M., Turnshek, D. A., & Wolfe, A. M. 1987, ApJ, 
322, 739 

Law, D. R., Steidel, C. C, Shapley, A. E., Nagy, S. R., Reddy, 

N. A., & Erb, D. K. 2011, ArXiv c-prints 
Lopez, G., & Chen, H.-W. 2011, ArXiv e-prints 
Lovegrove, E., & Simcoe, R. A. 2011, ApJ, 740, 30 
Marchesini, D., et al. 2007, ApJ, 656, 42 

Menard, B., Wild, V., Nestor, D., Quidcr, A., Zibctti, S., Rao, S., 

& Turnshek, D. 2011, MNRAS, 417, 801 
Mouhcine, M., Lewis, I., Jones, B., Lamareillc, F., Maddox, S. J., 

& Contini, T. 2005, MNRAS, 362, 1143 
Murdoch, H. S., Hunstead, R. W., Pettini, M., & Blades, J. C. 

1986, ApJ, 309, 19 
Nestor, D. B., Johnson, B. D., Wild, V., Menard, B., Turnshek, 

D. A., Rao, S., & Pettini, M. 2011, MNRAS, 412, 1559 
Nestor, D. B., Turnshek, D. A., & Rao, S. M. 2005, ApJ, 628, 637 
Noterdaeme, P., Srianand, R., & Mohan, V. 2010, MNRAS, 403, 

906 

Oppenheimer, B. D., Dave, R., Keres, D., Fardal, M., Katz, N., 
Kollmeier, J. A., & Weinberg, D. H. 2010, MNRAS, 406, 2325 

Papovich, C, Dickinson, M., Giavalisco, M., Conselice, C. J., & 
Ferguson, H. C. 2005, ApJ, 631, 101 

Poli, F., et al. 2003, The Astrophysical Journal Letters, 593, LI 

Prochaska, J. X., O'Meara, J. M., & Worseck, G. 2010, ApJ, 718, 
392 

Prochter, G. E., Prochaska, J. X., , & Buries, S. M. 2006, The 

Astrophysical Journal, 639, 766 
Quider, A. M., Nestor, D. B., Turnshek, D. A., Rao, S. M., 

Monier, E. M., Weyant, A. N., & Buschc, J. R. 2011, AJ, 141, 

137 

Roddy, N. A., Steidel, C. C, Pettini, M., Adelberger, K. L., 
Shapley, A. E., Erb, D. K., & Dickinson, M. 2008, ApJS, 175, 
48 



Ribaudo, J., Lehner, N., & Howk, J. C. 2011, ApJ, 736, 42 
Rubin, K. H. R., Weiner, B. J., Koo, D. C, Martin, C. L., 

Prochaska, J. X., Coil, A. L., & Newman, J. A. 2010, ApJ, 719, 

1503 

Sargent, W. L. W., Steidel, C. C, & Boksenberg, A. 1988, ApJ, 
334, 22 

Schechter, P. 1976, ApJ, 203, 297 
Schneider, D. P., et al. 2010, AJ, 139, 2360 
Simcoe, R. A. 2011, ApJ, 738, 159 

Simcoe, R. A., ct al. 2008, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Scries, Vol. 7014, 
Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Scries 

Simcoe, R. A., ct al. 2010, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 7735, 
Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Scries 

Simcoe, R. A., et al. 2011, ArXiv e-prints 

Steidel, C. C, Dickinson, M., & Pcrsson, S. E. 1994, ApJ, 437, 
L75 

Steidel, C. C, Kollmeier, J. A., Shapley, A. E., Churchill, C. W., 
Dickinson, M., & Pettini, M. 2002, ApJ, 570, 526 

Steidel, C. C, & Sargent, W. L. W. 1992, ApJS, 80, 1 

Tinker, J. L., & Chen, H. 2010, ApJ, 709, 1 

Tytler, D., Boksenberg, A., Sargent, W. L. W., Young, P., & 
Kunth, D. 1987, ApJS, 64, 667 

Vacca, W. D., Cushing, M. C, & Rayner, J. T. 2003, PASP, 115, 
389 

van Dokkum, P. G., et al. 2010, ApJ, 709, 1018 
Weiner, B. J., et al. 2009, ApJ, 692, 187 

Weymann, R. J., Williams, R. E., Peterson, B. M., & Turnshek, 

D. A. 1979, ApJ, 234, 33 
Willmer, C. N. A., et al. 2006, ApJ, 647, 853 
Wolf, C, Meisenheimer, K., Rix, H.-W., Borch, A., Dye, S., & 

Kleinheinrich, M. 2003, A&A, 401, 73 
Zhu, G., Moustakas, J., & Blanton, M. R. 2009, ApJ, 701, 86 
Zibetti, S., Menard, B., Nestor, D. B., Quidcr, A. M., Rao, S. M., 

& Turnshek, D. A. 2007, ApJ, 658, 161 



